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Abstract 

The conformal structure of the Schwarzschild-de Sitter spacetime is analysed using the 
extended conformal Einstein field equations. To this end, initial data for an asymptotic initial 
value problem for the Schwarzschild-de Sitter spacetime is obtained. This initial data allows 
to understand the singular behaviour of the conformal structure at the asymptotic points 
where the horizons of the Schwarzschild-de Sitter spacetime meet the conformal boundary. 

Using the insights gained from the analysis of the Schwarzschild-de Sitter spacetime in a 
conformal Gaussian gauge, we consider nonlinear perturbations close to the Schwarzschild- 
de Sitter spacetime in the asymptotic region. We show that small enough perturbations of 
asymptotic initial data for the Schwarzschild de-Sitter spacetime give rise to a solution to the 
Einstein held equations which exists to the future and has an asymptotic structure similar 
to that of the Schwarzschild-de Sitter spacetime. 
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1 Introduction 

The stability of black hole spacetimes is, arguably, one of the outstanding problems in mathemat¬ 
ical General Relativity. The challenge in analysing the stability of black hole spacetimes lies in 
both the mathematical problems as well as in the physical concepts to be grasped. By contrast, 
the nonlinear stability of Minkowski spacetime —see e.g. mm- and de Sitter spacetimes —see 
umiH]— are well understood. 

The results in mm show that the so-called conformal Einstein field equations are a pow¬ 
erful tool for the analysis of the stability and global properties of vacuum asymptotically simple 
spacetimes —see dUllllllllig. They provide a system of field equations for geometric objects 
defined on a 4-dimensional Lorentzian manifold (A4,g), the so-called unphysical spacetime^ which 
is conformally related to a spacetime (A4,g), the so-called physical spacetime, satisfying the Ein¬ 
stein field equations. The conformal framework allows to recast global problems in {M.,g) as 
local problems in (A4,g). The metrics g and g are related to each other via a rescaling of the 
form g = ’Efg where S is a so-called conformal factor. Crucially, the conformal Einstein field 
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equations are regular at the points where S = 0 —the so-called conformal boundary. Moreover, 
a solution thereof implies, wherever S ^ 0, a solution to the Einstein field equations. 

At its core, the conformal Einstein field equations constitute a system of differential conditions 
on the curvature tensors respect to the Levi-Civita connection of g and the conformal factor S. 
The original formulation of the equations as given in, say [TTl [13] , requires the introduction of so- 
called gauge source functions. An alternative approach to gauge fixing is to adapt the analysis to a 
congruence of curves. In the context of conformal methods, a natural candidate for a congruence 
is given by conformal geodesics —see pSl [53] . To combine gauges based on the properties of 
congruences of conformal geodesics with the conformal Einstein field equations, one needs a more 
general version of the latter —the so-called extended conformal Einstein field equations |2()j . The 
extended conformal field equations have been used to obtain an alternative proof of the semiglobal 
nonlinear stability of the Minkowski spacetime and of the global nonlinear stability of the de- 
Sitter spacetime —see [30]. In view of these results, a natural question is whether conformal 
methods can be used in the global analysis of spacetimes containing black holes. This article 
gives a first step in this direction by analysing certain aspects of the conformal structure of the 
Schwarzschild-de Sitter spacetime. 


1.1 The Schwarzschild-de Sitter spacetime 


The Schwarzschild-de Sitter spacetime is a spherically symmetric solution to the vacuum Einstein 
field equations with Cosmological constant. It depends on two parameters: the Cosmological 
constant A and the mass parameter m. The assumption of spherical symmetry almost completely 
singles out the Schwarzschild-de Sitter spacetime among the vacuum solutions to the Einstein field 
equations with de Sitter-like Cosmological constant. The other admissible solution is the so-called 
Nariai spacetime. This observation can be regarded as a generalisation of Birkhoff’s theorem —see 
m- For small values of the areal radius r, the solution behaves like the Schwarzschild spacetime 
and for large values its behaviour resembles that of the de Sitter spacetime. In the Schwarzschild- 
de Sitter spacetime the relation between the mass and Cosmological constant determines the 
location of the Cosmological and black hole horizons. 

The presence of a Cosmological constant makes the Schwarzschild-de Sitter solution a con¬ 
venient candidate for a global analysis by means of the extended conformal field equations: the 
solution is an example of a spacetime which admits a smooth conformal extension towards the 
future (respectively, the past) —see Figures and in the main text. This type of spacetimes 
are called future (respectively, past) asymptotically de Sitter —see Section 2.1 for definitions and 
HI [29] for a more extensive discussion. As the Cosmological constant takes a de Sitter-like value, 
the conformal boundary of the spacetime is spacelike and, moreover, there exists a conformal rep¬ 
resentation in which the induced 3-metric on the conformal boundary ..f is homogeneous. Thus, 
it is possible to integrate the extended conformal field equations along single conformal geodesics. 

In this article we analyse the Schwarzschild-de Sitter spacetime as a solution to the extended 
conformal Einstein field equations and use the insights thus obtained to discuss nonlinear pertur¬ 
bations of the spacetime. A natural starting point for this discussion is the analysis of conformal 
geodesic equations on the spacetime. The results of this analysis can, in turn, be used to rewrite 
the spacetime in the conformal gauge associated to these curves. However, despite the fact that 
the conformal geodesic equations for spherically symmetric spacetimes can be written in quadra¬ 
tures [23]) in general, the integrals involved cannot be solved analytically. In view of this difficulty, 
in this article we analyse the conformal properties of the exact Schwarzschild-de Sitter spacetime 
by means of an asymptotic initial value problem for the conformal field equations. Accordingly, 
we compute the initial data implied by the Schwarzschild-de Sitter spacetime on the conformal 
boundary and then use it to analyse the behaviour of the conformal evolution equations. An 
important property of these evolution equations is that their essential dynamics is governed by 
a core system. Consequently, an important aspect of our discussion consists of the analysis of 
the formation of singularities in the core system. This analysis is irrespective of the relation 
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Figure 1: Schematic depiction of the Main Result. Development of asymptotic initial data close 
to that of the Schwarzschild-de Sitter spacetime in the global representation —the initial metric 
is h, the standard metric on S^, and the asymptotic points Q and Q' are excluded (denoted by 
empty circles in the diagram). Figures a), b) and c) illustrate the evolution of initial data close 
to the Schwarzschild-de Sitter spacetime in the subextremal, extremal and hyperextremal cases 
respectively. See also Figures [3|4| and [5} 
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Figure 2: Schematic depiction of the Main Result. Development of asymptotic initial data close 
to that of the Schwarzschild-de Sitter spacetime in the representation in which Theorem was 
obtained. The initial metric is h, the standard metric on K x S^, and the asymptotic points Q and 
Q! are at infinity respect to h —since h and h are conformally flat one has h = oPh. The initial 
data for the subextremal, extremal and hyperextremal cases is formally identical. Small enough 
perturbations the development has the same asymptotic structure as the reference spacetime 


between A ^ 0 and m. This allows us to formulate a result which is valid for the subextremal, 
extremal and hyperextremal Schwarzschild-de Sitter spacetime characterised by the conditions 
0 < 9m^|A| < 1, 9m^|A| = 1 and 9 to^|A| > 1 respectively. 

1.2 The main result 

The analysis of the conformal properties of the Schwarzschild-de Sitter spacetime allows us to 
formulate a result concerning the existence of solutions to the asymptotic initial value problem 
for the Einstein field equations with de Sitter-like Cosmological constant which can be regarded 
as perturbations of the asymptotic region of the Schwarzschild-de Sitter spacetime —see Figures 
[Hand [21 Our existence result can be stated as: 

Main Result {asymptotically de Sitter spacetimes close to the asymptotic region of 
the SdS spacetime). Given asymptotic initial data which is suitably close to data for the 
Schwarzschild-de Sitter spaeetime there exists a solution to the Einstein field equations which exists 
towards the future (past) and has an asymptotic structure similar to that of the Schwarzschild-de 
Sitter spacetime —that is, the solution is future (past) asymptotically de Sitter. 

Remark 1. A detailed formulation of the Main Result of this article can be found in Section IT^ 
—see Theorem [H 

Our analysis of the conformal evolution equations governing the dynamics of the background 
solution provides explicit minimal existence intervals for the solutions. These intervals are cer¬ 
tainly not optimal. An interesting question related to the class of solutions to the Einstein field 
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equations obtained in this article is to obtain their maximal development. To address this prob¬ 
lem one requires different methods of the theory of partial differential equations and it will be 
discussed elsewhere. A schematic depiction of the Main Result is given in Figure 

As part of the analysis of the background solution we require asymptotic initial data for the 
Schwarzschild-de Sitter spacetime. The construction of this initial data allows to study in detail 
the singular behaviour of the conformal structure of the family of background spacetimes at the 
asymptotic points Q and Q', where the horizons of the spacetime meet the conformal boundary. 
As a consequence of the singular behaviour of the asymptotic initial data, the discussion of 
the asymptotic initial value problem has to exclude these points. In view of this, it turns out 
that a more convenient conformal representation to analyse the conformal evolution equations 
for both the exact Schwarzschild-de Sitter spacetime and its perturbations is one in which the 
the conformal boundary is metrically K x rather than S^\{Q, Q'} so that the problematic 
asymptotic points are sent to infinity —see Figure]^ In this representation, the methods of the 
theory of partial differential equations used to analyse the existence of solutions to the conformal 
evolution equations implicitly impose some decay conditions at infinity on the perturbed initial 
data. 

1.3 Related results 

The properties of the Schwarzschild-de Sitter spacetime have been systematically probed by means 
of an analysis of the solutions of the scalar wave equation using vector field methods —see [48j . 
This type of analysis requires special care when discussing the behaviour of the solution close to 
the horizons. In the asymptotic initial value problem considered in this article, the domain of 
influence of the initial data is contained in the region corresponding to the asymptotic region of 
the Schwarzschild-de Sitter spacetime. 

The properties of the Nariai spacetime —the other solution appearing in the generalisation 
of Birkhoff’s theorem to spacetimes with a de Sitter-like Cosmological constant— have been 
analysed by means of both analytic and numerical methods in BM- In particular, in the former 
reference it is shown that the Nariai solution does not admit a smooth conformal extension —see 
also |26| . Thus, it cannot be obtained from an asymptotic initial value problem. 

Finally, it is pointed out that the singularity of the Schwarzschild-de Sitter spacetime is 
not a conformal gauge singularity since CabcdC°‘^^‘^ —)■ oo as r —)• 0. Accordingly, theory of 
the extendibility of conformal gauge singularities as developed in [39] cannot be applied to our 
analysis. For any of the possible conformal gauges available, one either has a singularity of the 
Weyl tensor arising at a finite value of the parameter of a conformal geodesic or one has an 
inextendible conformal geodesic along which the Weyl tensor is always smooth. 

Notations and conventions 

The signature convention for (Lorentzian) spacetime metrics is (-I-, —, — , —In these conventions 
the Cosmological constant A of the de Sitter spacetime takes negative values. Cosmological 
constants with negative values will be said to be de Sitter-like. 

In what follows, the Latin indices a,b,c,--- are used as abstract tensor indices while the 
boldface Latin indices a, b^c, ■ ■ ■ are used as spacetime frame indices taking the values 0,..., 3. 
In this way, given a basis {ea} a generic tensor is denoted by Tab while its components in the given 
basis are denoted by Tab = Tabea°'eb^. We reserve the indices i, j, fc,... to denote frame spatial 
indices respect to an adapted frame taking the values 1, 2, 3. We make systematic use of spinors 
and follow the conventions and notation of Penrose & Rindler |45] — in particular, a, b, Ci ■ ■ ■ 
are abstract spinorial indices while a, b, Ct ■ ■ will denote frame spinorial indices with respect 
to some specified spin dyad {ca^}- Our conventions for the curvature tensors are fixed by the 
relation: 

(V,V, - Vb^aX = R^dabV’^- 
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In addition, D^{A), H{A), J^{A) and I^{A) will denote, respectively, the future (past) 
domain of dependence, the Cauchy horizon, the causal and the chronological futures (pasts) of A 
—see e.g. [Ml 155] . 


2 The asymptotic initial value problem in General Relativ¬ 
ity 

In this section we briefly revisit the notion of asymptotically de Sitter spacetimes —see (TJ [Ml |M] . 
After that, we review the properties of the extended conformal Einstein field equations that 
will be used in our analysis of the Schwarzschild-de Sitter spacetime. This general conformal 
representation of the Einstein field equations was originally introduced in |5D] — see also [MIIMIIM] 
for further discussion. For completeness, the conformal constraint equations are presented —see 
m m dzi ng. In addition, we provide a discussion on the notion of conformal geodesics and 
conformal Gaussian systems of coordinates —see (JHl HH IMl HO] • In this section we also discuss 
how to use the conformal field equations expressed in terms of a conformal Gaussian system to 
set up an asymptotic initial value problem for a spacetime with a spacelike conformal boundary. 
We conclude this section with a discussion of the structural properties of the conformal evolution 
equations in the framework of the theory of symmetric hyperbolic systems contained in |S7|. 

2.1 Asymptotically de Sitter spacetimes 

A spacetime (AI,g) satisfying the vacuum Einstein field equations 

^ab — ( 1 ) 

is future asymptotically de Sitter if there exist a spacetime with boundary (AI,g), a smooth 
conformal factor S and a diffeomorphism Lp ■. M. C AI, such that: 

S > 0 in U, 

S = 0 and dS 7 ^ 0 on = dU, 

is spacelike—i.e. g(dS, dS )>0 on , 

lies to the future of JA —i.e. C 

Observe that this definition does not restrict the topology of In particular, it does not have 
to be compact —see [29] . The notion of past asymptotically de Sitter is defined in analogous way. 
Additionally, (A4,g) is asymptotically de Sitter if it is future and past asymptotically de Sitter. 
Notice that a spacetime which is asymptotically de Sitter is not necessarily asymptotically simple 
—see |36j for a precise definition of asymptotically simple spacetime. In the following, in a slight 
abuse of notation, the mapping p : Ai ^ U Ai will be omitted in the notation and we write 

g = ( 2 ) 

Furthermore, the term asymptotic region will be used to refer to the set of a future 

asymptotically de Sitter spacetime or of a past asymptotically de Sitter spacetime. 

2.2 The extended conformal Einstein field equations 

In this section we provide a succinct discussion of the extended conformal Einstein field equations. 
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2.2.1 Basic notions 


Given any connection V over a spacetime manifold At, the torsion and Riemann curvature 
tensors are defined, respectively, by the expressions 

i'^a'^b — Vf,Va)(^ = Iia‘^b'^c4’: a'^ b — a)u^ = R'^ dabu'^ + , 

where 4> and u‘^ are smooth scalar and vector fields respectively, while Sa'^b and R'^cab denote the 
torsion and Riemann tensors of V. 


2.2.2 Frames and connection coefficients 

Let {ea} denote a set of frame fields on At and let {ti;“} be the associated coframe. One has 
that (a;“, 65 ) = Sb°'- We define the frame metric as gab = oi^a, 65 ) —in abstract index notation 
gab = Ga°‘eb^gab- From now on we will restrict our attention to orthonormal frames, so that 
gab = Vab, where consistent with our signature conventions gab = diag(l, — 1, — 1, — 1). The 
metric g is then expressed in terms of the coframe {uj°'} as 

g = gabba°' ® u>^. 

The connection coefficients Ta^^b of the connection V with respect to the frame {£„} are 
defined via the relation 

^a^b — Fq bCc, 

where Va = ea“Va denotes the covariant directional derivative in the direction of Ba- The torsion 
of V can be expressed in terms of the frame {ea} and the connection coefficients Fa'^b via 

Sa'^bGc = [ea,eb] - (F„^b - F5^„)ec. 


2.2.3 Conformal rescalings 


Following the notation introduced in Section 2.1 two spacetimes {M,g) are said to be {M,g) 


conformally related if the metrics g and g satisfy equation ([^ for some scalar field S. In the 
remainder of this article the symbols V and V will be reserved for the Levi-Civita connection of 
the metrics g and g. The connection coefficients of V and V are related to each other through 
the expression 


a b 


= F 


a b 


^ab 


cd'' 


where 

Sab‘^'^ = Sa^^Sb^ + Sb^^Sa^'- gabg'^'^ and T„ = S-iV„S. 
In particular, observe that the 1-form T = is exact. 


2.2.4 Weyl connections 

A Weyl connection V is a torsion-free connection satisfying the relation 

^ agbc — ‘^fagbc^ (3) 

where fa is an arbitrary 1 -form —thus, V is not necessarily a metric connection. Property ([^ 
is preserved under the conformal rescaling as it can be verified that Vattbc = —‘^faflbc where 
fa = fa + '^a- The Connection coefficients of V are related to those of V through the relation 

fa"b = r„"b + Sab^'^fd. (4) 
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A Weyl connection is a Levi-Civita connection of some element of the conformal class \g\ if and 
only if the 1-form is exact. The Schouten tensor Lab of the connection V is defined as 


Lab — i2^ah 

The Schouten tensors of the connections V and V are related to each other by 


1 


(5) 


Lab - Lab = ^afb " ^Sab^^fcfd 
Notice that, in general, Lab ^ L(^ab)- 

2.2.5 The extended conformal Einstein field equations 

From now on, we will consider Weyl connections V related to a conformal metric g as in equation 
(§. Let P°‘bcd denote the geometric curvature of V —that is, the expression of the Riemann 
tensor of V written in terms of derivatives of the connection coefficients Ta'^f,- 

P°'bcd = GaCfb^d) — ebita^^d) + Lf^^diLb^a — L b) + L b^ dL f — L d^ b^^ f ■ 

The expression of the irreducible decomposition of Riemann tensor R°'bcd given by 

p“bcd = bed + ‘^Sb[c°'^ Lb]f (6) 

will be called the algebraic curvature. In the last expression dPbed represents the so-called rescaled 
Weyl tensor, defined as 

ja _ -r —1/^a 

bcA — 


bcdi 


where C°’bcd is the Weyl tensor of the metric g. Despite the fact that the definition of the rescaled 
Weyl tensor may look singular at the conformal boundary, it can be shown that under suitable 
assumptions the tensor d°‘bcd it is regular even when S = 0. Finally, let us introduce a 1-form d 


defined by the relation 

— ^fa H" ^ 


With the above definitions one can write the vacuum extended conformal Einstein field 

equa- 

tions as 



Ea"b 

= 0, ^“bed = 0, Acdb = 0, Abed = 0 

(7) 

where 

Ea'^b = [eajGb] — (f'a'^b — fb‘^a)ec. 

(8a) 


_ pa pa 

“ bed — ^ bed bed-) 

(8b) 


^edb — ^ e^bd d^eb dad bed — 0: 

(8c) 


^bed — ^ad bed fad bed- 

(8d) 


The fields 


bed 


, A and Abed will be called zero-quantities. 


The geometric meaning of the extended conformal field equations is the following: fia^b = 0 
describes the fact that the connection V is torsion-free. The equation E°'bcd = 0 expresses the 
fact that the algebraic and geometric curvature coincide. The equations Acdb = 0 and Abed = 0 
encode the contracted second Bianchi identity. Observe that there is no differential condition for 
neither the 1-form d nor the conformal factor. In Section 12.2.61 it will be discussed how to fix 
these fields by gauge conditions. 




In order to relate the conformal equations 0 to the vacuum Einstein held equations ( |41| ) one 
introduces the constraints 

(5a = 0, 7ab = 0, Cab = 0 (9) 

encoded in the supplementary zero-quantities 


Sa = da- Efa - VS, 

7ab = 0^^"^ Vab Va(‘^ db) 


’ 3 ' —2 c* cd 

*^1 


dc^dt 


Cab — d^[ab] ^ a fb- 

The hrst equation in ([^ encodes the dehnition of the 1-form da', the second equation in ^ arises 
from the transformation law between the Schouten tensor Lab of V and the physical Schouten 
tensor Lab — \ f]ab determined by the Einstein held equations (411; the last equation in ([^ relates 
the antisymmetry of the Schouten tensor Lab to the derivative of the 1-form fa- 

The precise relation between the extended conformal Einstein held equations and the Einstein 
held equations is given by the following lemma: 

Lemma 1. Let Ta^c, ^ab, d“bcd) denote a solution to the vacuum extended conformal Ein¬ 
stein field equations ([^ for some choice of gauge fields CE^da) satisfying the constraint equations 
0 - Assume, further, that 


^0 


and 


det{'q°'^eaeb) 7 ^ 0 


on an open subset U C A4. Then 


9 = 0 ^Vab<^°' ® 


where {w“} is the coframe dual to {ea}, is a solution to the vacuum Einstein field equations (41) 
on lA. 

The proof of this lemma can be found in [MIISS]- 


2 . 2.6 Conformal geodesics and conformal Gaussian systems 

A conformal geodesic on a spacetime (M,g) consists of a pair {x{t), (3{t)) where x{t) is a curve 
with tangent x(t) and /3 (t) is a 1-form dehned along x{t) satisfying the conformal geodesic 
equations 


= -x‘^x^Sdb‘'^ fif, 
x^^cfia = -\x^Sj<^l3bfid + LaaX^ 


(11a) 

(lib) 


where Lr„ denotes the Schouten tensor of V and 


cd 


Sab^’^ = Sa^db'^ + da^db'^ - gabq 

In addition, it is convenient to consider a Weyl propagated frame —that is, a frame held {£„“} 
satisfying 

= -SacfiUa'^X^fif. 

The dehnition of conformal geodesics is motivated by the transformation laws of equations 
(11a l-( 11b) under conformal rescalings and transitions to Weyl connections. More precisely, given 
an arbitrary 1-form / one can construct a Weyl connection V as in equation Then, dehning 
(3 = (3 — f the pair {x{t), I 3 {t)) will satisfy the equations 


= -x'^x'‘Sdb°‘^Pf, 
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i"Vc/3a = 

where Lea is the Schouten tensor of the connection V as defined in equation ([^. If one chooses 
a Weyl connection V whose defining 1-form f coincides with the 1-form f3 of the V-conformal 


geodesic equations (lla)-(llb), then the conformal geodesic equations reduce to 

= 0 , 


Labi'' = 0 . 


( 12 ) 


Similarly, the Weyl propagation of the frame becomes 


= 0 . 


(13) 


The conformal geodesics equations admit more general reparametrisations than the usual 
affine parametrisation of metric geodesics. This is summarised in the following lemma: 

Lemma 2. The admissible reparametrisations mapping (non-null) conformal geodesics into (non¬ 
null) conformal geodesics are given by fractional transformations of the form 

ar -\-b 

T , 

cr -I- a 


where a, b,c,dG K. 


The proof of this lemma can be found in [21] — see also [521 ESj ■ Conformal geodesics allow 
to single out a canonical representative of the conformal class [g]. This observation is contained 
in the following key result: 


Lemma 3. Let {A4,g) be a spacetime where g is a solution to the vacuum Einstein field equations 
(41). Moreover, let [x(t), (3(t)) satisfy the conformal geodesic equations (llal-(llb) and let {e a} 
denote a Weyl propagated g-orthonormal frame along x(t) with 


a = 


such that 

g{x,x) = 1. 

Then the conformal factor 0 is given, along x{t), by 

0(r) = 0*-b 0*(r - T*) -h ^0 *(t - T*)^ (14) 

where the coefficients 0* = 0(r*), 0* = 0(t*) and = 0(t*) are constant along the conformal 
geodesic and satisfy the constraints 

0* = (/3*,a;,^)0*, 0*0* = ^^“(/J*,/?*)-b ^A. (15) 

Moreover, along each conformal geodesic 

0/3o = 0, 0/3i = 0*/3i*, 


where Pa = {(^,ea)- 

The proof of this Lemma and a further discussion of the properties of conformal geodesics can 
be found in [201ES] . 

For spacetimes with a spacelike conformal boundary the relation between metric geodesics 
and conformal geodesics is particularly simple. This observation is the content of the following: 
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Lemma 4. Any conformal geodesic leaving ) orthogonally into the past (future) is up 

to reparametrisation a timelike future (past) complete geodesic for the physical metric g. The 
reparametrisation required is determined by 


df 

dr 


1 


0(t) 


(16) 


where f is the g-proper time and t is the g-proper time and g = Q^g. 
The proof of this Lemma can be found in [55]. 


2.2.7 Conformal Gaussian systems 

In what follows it will be assumed that there is a region of the spacetime (At,g) which can 
be covered by non-intersecting conformal geodesics emanating orthogonally from some initial 
hypersurface S. Using Lemma the conformal factor (141 is a priori known and completely 
determined from the specification of 0*, 0* and 0* on S. A conformal Gaussian system is then 
constructed by adapting the time leg of the g-orthonormal tetrad {£„} to the tangent to the 
conformal geodesic {x{t) , ( 3 {t)) —i.e. one sets eo = ds. The rest of the tetrad is then assumed to 
be Weyl propagated along the conformal geodesic. If one writes this condition together with the 
conformal geodesic equations expressed in terms of the Weyl connection singled out by ( 3 , as in 


equations (12) and (13), one obtains the gauge conditions 


0 b 


= 0, 


Loa — 0 , 


/o =0. 


(17) 


One can further specialise the gauge by using the parameter r along the conformal geodesics as 
a time coordinate so that 

eo = dr- (18) 

Now, consider a system of coordinates (r, x“) where (x“) are some local coordinates on S. 
The coordinates (a;“) are extended off the initial hypersurface S by requiring them to remain 
constant along the conformal geodesic which intersects a point p G S with coordinates (a;“). This 
type of coordinates will be called a conformal Gaussian coordinate system. This construction 
naturally leads to consider a 1-1-3 decomposition of the field equations. 


2.2.8 Spinorial extended conformal Einstein field equations 


A spinorial version of the extended conformal Einstein held equations (8a)-(8d) is readily obtained 


by suitable contraction with the Infeld-van der Waerden symbols (7°' aa' ■ Given the components 
Tab^ of a tensor the components of its spinorial counterpart are given by 


Taa' 


BB' 


CC 


= Tab^ <7 AA' 


bOC 

<7bB' (7 c 


where, 


0 — 

<7 A A' = 


<7 A A' = 


^AA! 

(7 0 


AA' 
a 2 


72 VO 

-(" 

72 Vi 

— r 

■ 72 VO 

i 

72 V-i 0 


0 

—1 

0 

0 

1 


7AA' 


(7 A A' = 


^AA! 
(7 1 


„AA' 

<7 3 


^(0 

72 VI 

72 VO 
^j_/o 
“ 72 VI 

= -0 

“ 72 VO 


1 

0 

0 

-1 

1 

0 

0 

-1 
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In particular, the spinorial counterpart of the frame metric gab = Vab is given by gAA'BB' = 
cab^A'B' while the frame {e^} and coframe imply a frame {baa'} and a coframe 
such that 

gi^AA' ,eBB') = ^AB^A'B'- 


If one denotes with the same kernel letter the unknowns of the frame version of the extended 
conformal Einstein field equations one is lead to consider the following spinorial zero-quantities: 

^AA'BB' = \eAA’.^BB’] — AA' BB’—^BB’ AA')^CC'■, (19a) 

•hCC — t-CC hCC 

^ DD'AA'BB' = r DD'AA'BB' — B DD'AA'BB', (190) 

^CC'DD'BB' = y CC'LdD'BB' — ^DD'LcC'BB' — dAA'd^^ BB'DD', (19c) 
^BB'CC'DD' = AA'd^'^ BB'CC'DD' — fAA'd'^^ BB'CC'DD' ■ (19d) 


The spinorial version of the extended conformal Einstein field equations are then succinctly 
written as 

^AA'BB' = 0, 5*"^ DD'AA'BB' = 0, ^CC'DD'BB' = 0, ^BB'CC'DD' = 0. (20) 


In the spinor description one can exploit the symmetries of the fields and equations to obtain 
expressions in terms of lower valence spinors. In particular, one has the decompositions 


dAA'BB'CC'DB' = —(t>ABCD^A'B'^C'D' — (f’A'B'Cn'iABiCB, 

•p BB' p S _ B' I p B' _ B 

^ AA' CC' = ^ AA' C^C +^AA' C'^C , 

where 4>abcd = 4>(Abcd) are the components of the rescaled Weyl spinor and TaA'^ c are the 
reduced connection coefficients of V. Using the spinorial version of equation Q and contracting 
appropriately one obtains 

TcC'AB = TcC'AB — ^AcfBC- (21) 

Likewise, one has the following reduced curvature spinors 


B dAA'BB 


, = -R^^' 
2 


DQ'AA'BB' 


= ^AA’ — esB' (J'aA'^L>) 

f Cf F f Cf F',f Cf F 
— t FB' Dt aA' B — t BF' Dt AA' B' + t FA' BB' A 

I p Op F' iP ^ P ^P ^ 

-rl AF' BB' A' -r ^ AA' BB' D — ^ BB' AA' D, 


Pabcc'dd' = 2^^^ 


BQ'CC'DD' 


= —GXjjABCDeC'D' — Lc'{AB)D'^CE> — 2^AB {LcC'DB' — LdD'CC' 


^AA'BC = 2^AA'BQ'C 


Q' 


= a4>ABCQ - A'4>ABCQ- 


With these definitions, the spinorial extended conformal Einstein field equations can be alterna¬ 
tively written as 

^AA'BB' = 0, DA A'BB' = 0, ^CC'DD'BB' = 0, ^BB'CD = 0. (22) 


The last set of equations is completely equivalent to the equations in (20). 
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2.2.9 Space spinor formalism 

In what follows, let the Hermitian spinor denote the spinor counterpart of the vector •\/2eo“. 
In addition, let {ea^} with ei^ = denote a spinor dyad such that 


^AA' ^ A^ A' 
T = Co eo' 


Cl ei' ■ 


We have chosen the normalisation taa' = 2, in accordance with the conventions of [IH]- In 
what follows let denote the components of respect to {c^a}- The Hermitian spinor 

can be used to perform a space spinor split of the frame {caa'} and coframe }. 

Namely, one can write 


^ A A' Ft 

e-AA'= 7,'^ e-T a'Sab, 


, A' nA 

UJ ~ 2^ U) + Tc w , 


(23) 


where 


_ pp' _ p' _ pp' 

e = T epp', eAB=T(^A ^b)P'j u3 = tpp'U3 


, ,AB (A , ,B)P' 

LJ = —T' piUJ ' 


It follows from the above expressions that the metric g admits the split 

g = (g) o) + Kabcd^^^ ® 


where 


hABCD = g{eAB,eCL>) = —^A(C^D)B- 


Similarly, any general connection V can be split as 


Vaa' = ■^taa''P — ta'^'Daq, 

where 

V = V aa' and T^ab = T(^b^a}A', 

denote, respectively, the derivative along the direction given by and Vab is the Sen con¬ 

nection of V relative to . 

The Hermitian spinor induces a notion of Hermitian conjugation: given an arbitrary 

spinor with components ^aa' its Hermitian conjugate has components 

i4jd = tu^JTXa^ = (24) 


where the bar denotes complex conjugation. In a similar manner, one can extend the definition 
to contravariant indices and higher valence spinors by requiring that (tt/x)^ = 


2.3 Conformal evolution and constraint equations 

In this section the evolution equations implied by the extended conformal field equations and the 
conformal Gaussian gauge are discussed. In addition, a brief overview of the conformal constraint 
equations is given. 


2.3.1 Conformal Gauss gauge in spinorial form and evolution equations 

The space spinor formalism leads to a systematic split of the extended conformal Einstein field 
equations (22) into evolution and constraint equations. To this end, one performs a space spinor 
split for the helds baa', /aa', Taa', Taa'^c ■ 
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The frame coefficients eAA'°' satisfy formally identical splits to those in (23), where baa' = 
eAA'°'Ca with Ca G {dr,Ci} representing a fixed frame field —the latter is not necessarily g- 
orthonormal. Observe that, in terms of tensor frame components, the gauge condition (18) 
implies that eo“ = <5o“. The gauge conditions 0 and (|18[) are rewritten as 


AA.'B n 
^ j- AA' C = U, 


e A A' = V^dr, 

In addition, we define 

^ABCD = ^AA 'CD, ^ABCn = Tb^ ^AA'CD^ 

Labcd = Tb^ Tb^ LaA'CCG 


T-^-^ Laa'bb' = 0. 


= 'TB^^AA'CDi fAB = TB^fAA', 
^ABCD = Lab{CD) ^AB = LabQ^- 


(25) 


(26a) 

(26b) 


Recalling equation (21) one obtains 


^ABCD — ^ ABCD — ^CAfnA'TB 


A' 


where Tabcb = tb"^ ^AA'cn- This relation allows us to write the equations in terms of the 
reduced connection coefficients of the Levi-Civita connection of g instead of the reduced connec¬ 
tion coefficients of V. Only the spinorial counterpart of the Schouten tensor of the connection V 
will not be written in terms of its Levi-Civita counterpart. Exploiting the notion of Hermitian 
conjugation given in equation (p4|) one defines 


XABCD 




+ r 


t 

ABCn )> 


^ABCD = ^ {yABCD — T 


t 

ABCDI > 


where xabcd and ^abcd correspond, respectively, to the real and imaginary part of the connec¬ 
tion coefficients Tabcd- We define the electric and magnetic parts of the rescaled Weyl spinor 
as 


VABcn = X {(jiABcn + <I^ABcn) > ^-abcd = (<^ 


ABCD 


^^ABCD 


The gauge conditions (25) can be rewritten in terms of the spinors defined in (26a) as 
fAB = f{AB), = 0, Lq^aB=0- 

The last condition implies the decomposition 


(27) 


Labcd = Qabcd + ^scd&ab 

for the components of the spinorial counterpart of the Schouten tensor of the Weyl connection 
where Qabcd = L(^ab){cd) and Oab = Labq^- 

The fields defined in the previous paragraphs allow us to derive from the expressions 


T^^'tAA'^^' 


BB'Bpp'°' = 0, 
^AA'BB'CC = 0, 


^CC'G n 

T ^ABCC'DD' = U, 

^|A'|BCD) = 0, 


'BB'CC 

a set of evolution equations for the fields 

XABCD, ^ABCD, Bab° , Bab^ , JaB, QaBCD, ^AB, i’ABCD- 

Explicitly, one has that 

drS-AB^ = —X{AB)^‘^epQ° — /aB, 
dr^AB^ = —X{AB)^^^Pq'', 


(28a) 

(28b) 


(29a) 

(29b) 
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dr^ABCD = -X(AB)^‘^iPQCn + -^{<^ABX{BD)PQ + <^BDX(AC)PQ)f^'^, (29c) 

— '^XAB(C^fB)E — 2 {^AcQbD + £BB&Ac) — ^Qf^ABCD, (29d) 

drfAB = —X{AB)^'^fpQ + -^^AB, (29e) 

dTX(AB)CD = —XAB^^XPQCD — Qabcb + QiJABCD, (29f) 

dr^ABCD = —X(AB)^‘^ LpQi^CD) — Qvabcl) + (AfJ'B)CnP, (29g) 

dr^AB = —X(AB)^^^EF + V2d^^r]ABPQ, (29h) 

dT(t>ABCD — y^T^{A'^4>BCD)Q = 0. (29i) 


The following proposition relates the discussion of the conformal evolution equations and the 
full set of extended conformal field equations given by Q: 

Lemma 5 {propagation of the constraints and subsidiary system). Assume that the evo¬ 


lution equations extracted from equations (28a I-(28b) and the conformal Gauss gauge eonditions 


(27) hold. Then, the independent components of the zero quantities 


V BB' 
^AA' CC 


-ABCC'DEI' 


^AA'BB'CC 


^AA'i lAA'BB' 


Caa', 


which are not determined by the evolution equations or the gauge conditions, satisfy a symmetric 
hyperbolic system of equations whose lower order terms are homogeneous in the zero-quantities. 

The proof of Lemma can be found in uniiiiiissi — see also mi for a discussion of these 
equations in the presence of an electromagnetic field. 

The most important consequence of Lemma ^ is that if the zero-quantities vanish at some 


initial hypersurface and the evolution equations (|29a|)- (29h) are satisfied, then the full extended 


conformal Einstein field equations encoded in (20) are satished in the development of the ini¬ 


tial data. This is a consequence of the standard uniqueness result for homogeneous symmetric 
hyperbolic systems. 


2.3.2 Controlling the gauge 


The derivation of the conformal evolution equations (29a)-(29h) is based on the assumption of 


the existence of a non-intersecting congruence of conformal geodesics. To verify this assumption 
one has to analyse the deviation vector of the congruence. 


Let 2: denote the deviation vector of the congruence. One has then that 

[a;, 2:1 = 0. 


(30) 


Now, let 2;^"^ denote the spinorial counterpart of the components 2“ of 2; respect to a Weyl 
propagated frame {e^}. Following the spirit of the space spinor formalism one dehnes zab = 
tb^ zaa'- This spinor can be decomposed as 


1 

ZAB = 2^^AB 


■ Z(AB) 


The evolution equations for the deviation vector can be readily deduced from the commutator 


(30). Expressing the latter in terms of the fields appearing in the extended conformal field 


equations one obtains 


drZ = fABZ 


{AB) 


(31a) 
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dTZ(AB) = XCD(AB)Z^^^'^ 


(31b) 


The congruence of conformal geodesics is non-intersecting as long as z^ab) 7^ 0. Once one has 
solved equations (29a I- (29i) one can substitute Jab and xabcd into equation (311 and analyse 
the evolution of the deviation vector —for further discussion see HD] • 


2.3.3 The conformal constraint equations 


The conformal constraint equations encode the set of restrictions induced by the zero-quantities 
on the various fields on hypersurfaces of the unphysical spacetime {Ai,g). In what follows, we 
will consider a setting where the 1-form f vanishes. Accordingly, the initial data for the extended 
conformal evolution equations (29a)-( 29h ) and those for the standard conformal Einstein field 
equations are the same —see Appendix [pj Now, let S denote a 3-dimensional submanifold of M. 
The metric g induces a 3-dimensional metric h = (p*g on S, where : 5 —>■ is an embedding. 

Similarly, one can consider a 3-dimensional submanifold 5 of with induced metric h = 
such that 

h = 


where denotes the restriction of the conformal factor to the initial hypersurface S —in Section 
|2.2.6| this restriction is denoted by 0*. 

Let Ua and fia with = 0ha be, respectively, the g-unit and g-unit normals, so that 
n“na = fL°"na = 1 —in accordance with our signature conventions for an spacelike hypersurface. 
With these definitions, the second fundamental forms Xab = ha’^VcUb and Xab = ha^Vcfib are 
related by the formula 

Xab ~ ^{Xab T ^^ab) 

where S = n^'Va^- 

The conformal constraint equations are conveniently expressed in terms of a frame {ej} 
adapted to the hypersurface S —that is, the vectors span TS and, thus, are orthogonal 
to its normal. All the fields appearing in the constraint equations are expressed in terms of this 
frame. The conformal constraint equations are then given by: 


DiDjn = —YjXij — ^Lij + shij, (32a) 

DiE = Xi^Dk^ - flLi, (32b) 

DiS = —LiT, — XlLi, (32c) 

^iLjk — Edijfc -j- dlj^ij D XI (^Xik^j Xjkd^i) ^ (32d) 

DiLj — DjLi = diijD^Xl -\- Xi'^Ljk ~ Xj^^ik^ (32e) 

D’^dkij = x’^idjk - x'^jdik, (32f) 

D^dij = x^^dijk, (32g) 

d^jXki dDkXji — Xldijf^ -f hijLf^ (32h) 

lij = Xldij -|- Lij Xk ixij ^xdij) + XkiXj ^XkiX ) (32i) 

A = - 3E2 - SDkXlD'^Xl, (32j) 


where D is the Levi-Civita connection on {S,h), lij is the associated Schouten tensor, dijk = 
diojk, dij = diojo, Li = Loi and s is a scalar field on S —see Appendixfor the definition of s 
in context of the conformal Einstein field equations. 


2.3.4 Constraints at the conformal boundary 

The conformal constraint equations simplify considerably on hypersurfaces for which = 0. If 
this is the case then equations (32a)-(32i) reduce to 


shij = Ey. 


ij I 


(33a) 
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Di^ = 0, (33b) 

DiS = —LiS, (33c) 

^ (33ci) 

DiLj DjLi = \i ^jk Xj ^iki (33e) 

D^dkij = x'^idjk - x'^jdik, (33f) 

A = -3E2, (33g) 

D^dij = x^'^dijk, (33h) 

ddjXki dDkXji (33i) 

lij — Lij ~ xixij ~ ^xdij) + XkiXj ~ ~^XklX ^ij . (33j) 


A procedure for obtaining a solution for these equations has been given in [HI HO]- Direct 
algebraic manipulations yield 


S = 



s = Ek, 


2"' 


Xij — ; dji — DiK,^ 

dijk = -^-^Vijk, 


(34a) 

(34b) 


where k is an smooth scalar function on the initial hypersurface and yijk denotes the components 
of the Cotton tensor of the metric h. The only differential condition that has to be solved to 
obtain a full solution to the conformal constraint equations is 


D^dij = 0 , 


(35) 


where dij is a symmetric tracefree tensor encoding the initial data for the electric part of the 
rescaled Weyl tensor. 


2.4 The formulation of an asymptotic initial value problem 

In this section we show how the conformal Gaussian gauge can be used to formulate an asymptotic 
initial value problem for the extended conformal Einstein field equations. Thus, in the sequel we 
consider an initial hypersurface on which the conformal factor vanishes so that it corresponds to 
the conformal boundary of an hypothetical spacetime. Accordingly, this initial hypersurface will 
be denoted by y. 


2.4.1 The conformal boundary 


Following Lemma we can set, without lost of generality, t* = 0 on Moreover, it will be 
assumed that fa vanishes initially. Accordingly, we have the initial condition /3* = ©T^d©*. 
Recalling that d — Q(3, and and using the constraints in (15) of Lemma it readily 

follows, for the asymptotically problem (in which ©* = 0), that 



Moreover, using again that d = ©/3 and requiring A* to be orthogonal to Jd (so that a;* = Bq), 
we have do* = ©*• It follows that 
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The coefhcient 0* is fixed by the requirement s = Sk on —see [4j. From the definition of 

s and = Vn© it follows that 


.. = I jv 


.v-e + Ffie) = j(e„E“). + 


= ^??“''(eae50)* + ^0*(r„“o). 


(36) 


Taking into account that 0 and Sj vanish at we have that 77 “^(eQ,e 50 )* = 0*. Using 
the solution to the constraints given in (34al-( |34b| ) and exploiting the properties of the adapted 
orthonormal frame we have (ra“o)* = (Tj^o)* = (Xi*)* = = 3 k. Substituting into (36) and 

using that s* = 0 *k one gets 

0* = 0*K. 


Summarising, for an asymptotic initial value problem the conformal factor implied by the 
conformal Gaussian gauge is given by 


0 (r) = 



(37) 


The conformal factor given by equation (37) is, in a certain sense, Universal. It does not encode 
any information about the particular details of the spacetime to be evolved from ^. As such, 
it can be used to analyse any spacetime with de Sitter-like Cosmological constant as long as 
the spacetime has at least one component of the conformal boundary. If k 7 ^ 0 the conformal 
boundary has two components located at 


T = 0 


and r 


2 

K 


The first zero corresponds to the initial hypersurface J^. The physical spacetime corresponds 
to the region where 0 7 F 0. Therefore, the roots of 0 render two different regions of {M,g) 
corresponding to two different conformal representation of {Ai,g). One of these representations 
corresponds to the region covered by the conformal geodesics with t € [— 2 /|k|, 0 ] or r S [ 0 , 2 /|k|] 
and other corresponds to the region covered by the conformal geodesics with r S [ 0 , 00 ) or 
T G (— 00 , 0] depending on the sign of k. 


Remark 2. The discussion of the previous paragraphs is formal: the component of the conformal 
boundary given by t = — 2/k may not be realised in a specific spacetime. This is, in particular, 
the case of the extremal and hyperextremal Schwarzschild-de Sitter spacetimes in which the 
singularity precludes reaching the second conformal infinity —see Figure]^ 


2.4.2 Exploiting the conformal gauge freedom 

The conformal freedom of the setting allows us to further simplify the solution to the conformal 
constraint equations at J^. Given a solution to the conformal Einstein field equations associated 
to a metric g, it follows from the conformal covariance of the equations and fields that the 
conformally related metric g' = 'd'^g for some 1 } is also a solution. On an initial hypersurface S 
the latter implies implies h' = ‘d^h. From the definition of the field s —see Appendix and 
the conformal transformation rule for the Ricci scalar one has that 


s; = C^s*+C"(Vc^^)*(V"0)*. 

Thus, the condition s' = 0 can be solved locally for i?*. Accordingly, one chooses so that 
K = 0. In this gauge Xij and vanish and at In addition, the conformal factor 

reduces to 
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In this representation 0 has only one zero and the second component of the conformal boundary 
(if any) is located at an infinite distance with respect to the parameter r. 


2.5 The general structure of the conformal evolution equations 


One of the advantages of the hyperbolic reduction of the extended conformal Einstein field equa¬ 
tions by means of conformal Gaussian systems is that it provides a priori knowledge of the 
location of the conformal boundary of the solutions to the conformal field equations. Following 
the discussion in Section 2.2.7 the conformal geodesics fix the gauge through equations 0 and 
(18). The last condition corresponds to the requirement on the spacetime to possess a congru¬ 
ence of conformal geodesics and a Weyl propagated frame —i.e. equations jg and ( |l^ are 
satisfied. As already mentioned, the system of evolution equations (29a)- (29h) constitutes a 


symmetric hyperbolic system. This is the key property for analysing the existence and stability 
of perturbations of suitable spacetimes using the extended conformal Einstein field equations. 

To discuss the structure of the conformal evolution system in more detail, let e denote the 
components of the frame bab, T the independent components of xabcd and ^abcl>, and 4> 
the independent components of the rescaled Weyl spinor (j)ABCD- Then the evolution equations 


(29a)-(29h) can be written as 


drV = Kw -f Q(r)r^ -I- L(a;)</>, 

(I + A°{e))drcj) + A^dicf} = B(r), 


(38a) 

(38b) 


where v represents the independent components of the spinors in the conformal evolution equa¬ 
tions except for the rescaled Weyl spinor whose components are represented by cj). In addition, I 
is the 5x5 identity matrix, K is a constant matrix, Q, A*^, A^, and B are smooth matrix valued 
functions of its arguments and L(a:) is a matrix valued function depending on the coordinates. 
To have an even more compact notation let u = (v^cp). Consistent with this notation, let u 
denote a solution to the evolution equations (38a)-(38b) arising from data u* prescribed on an 


hypersurface S. The solution u will be regarded as the reference solution. Consider a general 
perturbation succinctly written as u = u -|- u. Equivalently, one considers 


e = e -I- e, 


r = r + r. 


cp = (p + cp. 


(39) 


Recalling that u is a solution to the conformal evolution equations (38a)-(38b) and making use 


of the split (39) one obtains that 


drV = Kii -I- Q(f -I- f )u -1- Q(f )'u -I- L{x)^, 

(I -I- A°(e -I- e))dr^ + (I + A°(e -|- e))dr^ + A*(e -f e)di(p + 
A*(e -f e)di^ = B(f -f f )^ + B(f + f)^. 


(40a) 


(40b) 


Equations (40a) and (40b) are read as equations for the components of the perturbed fields v and 


(p. These equations are in a form where the theory of first order symmetric hyperbolic systems in 
m can be applied to obtain a existence and stability result for small perturbations of the initial 
data u*. This requires however, the introduction of the appropriate norms measuring size of 
the perturbed initial data ii*. This general discussion will not be developed further, instead, we 
particularise this discussion in Section [4.3| introducing the appropriate norms required to analyse 
the Schwarzschild-de Sitter spacetime as an asymptotic initial value problem. 


3 The Schwarzschild-de Sitter spacetime and its conformal 
structure 

In this section we briefly review general properties of the Schwarzschild-de Sitter spacetime that 
will be relevant for the main analysis of this article. 
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3.1 The Schwarzschild-de Sitter spacetime 

The Schwarzschild-de Sitter spacetime is the spherically symmetric solution to the Einstein field 
equations 

Rah — ^9ab (It) 

with, in the signature conventions of this article, a negative Cosmological constant given in static 
coordinates {t, r, 0 , ip) by 

QSdS = F{r)dt (g) dt — F(r)“^dr (g) dr — r^cr, (42) 

where the function F{r) is given by 

F(r) = 1 - — + ^Ar^, (43) 

r 3 

and (T is the standard metric on the 2-sphere 

(T = d0 (g d0 -I- sin^ Odp (g di^, 


with t e (— 00 , 00 ), r e (0,oo), 9 G [0,7r], p G [0, 27r). This solution reduces to the de Sitter 
spacetime when m = 0 and to the Schwarzschild solution when A = 0. 

Remark 3. In the following, we will only consider the case m > 0 and we will always assume a 
de Sitter-like value for the cosmological constant A. 

The location of the roots of the polynomial r — 2m -I- |Ar^ are determined by the relation 
between m and A; whenever 0 < 9m^|A| < 1 this polynomial has two distinct positive roots rb,rc 
and a negative root r_ located at 


n = 


2 fa 

— z= cos — -I- 

yiAi Vs 


(5) 


—^ COS , 

V3 

2 fa 

— 1 = cos — -I- 
^ V3 


47r 

Y 


27r 

T 


where cos a = The positive roots 0 < rt < r^ correspond, respectively, to a black 

hole-like horizon and a Cosmological-like horizon. One can classify this 2-parameter family of 
solutions to the Einstein field equations depending on the relation between the parameters m and 
A . The subextremal Schwarzschild-de Sitter spacetime arises when the relation between m and 
A satisfies 

0 < Qm^lAI < 1. (44) 


If condition (44) holds, one can verify that F{r) > 0 for r;, < r < Tc while F(r) < 0 in the regions 
0 < r < rb and r > rc- Consequently, the solution is static for rb < r < rc —see [?]• The extremal 
Schwarzschild-de Sitter spacetime is obtained by setting 


|A| = l/9m^. 


(45) 


If the extremal condition (45) holds, then the black hole and Cosmological horizons degenerate 


into a single Killing horizon at r = 3m. Moreover, one has that F{r) < 0 for 0 < r < 00 so that 
the hypersurfaces of constant coordinate r are spacelike while those of constant t are timelike and 
there are no static regions. In the extremal case the function F{r) can be factorised as 


F{r) = - 


(r — 3m)^(r -|- 6m) 
27m^r 


(46) 
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Figure 3: Penrose diagram for the subextremal Schwarzschild-de Sitter spacetime. The excluded 
points Q, Q' represent asymptotic regions where the Cosmological horizon appear to meet . 


As discussed in Section 3.1 this region of the spacetime does not belong to . 


In the hyperextremal Schwarzschild-de Sitter spacetime one considers 

Qm^lAI > 1. (47) 

In this case one has again F{r) < 0 for 0 < r < oo so that similar remarks as those for the 
extremal case hold. The crucial difference with the extremal case is that in the hyperextremal 
case there are no horizons. Finally, at r = 0 it can be verified that the spacetime has a curvature 
singularity irrespective of the relation between m and A —in particular, the scalar CabedC°'^‘^‘^i 
with C‘^bed the Weyl tensor of the metric gsdS, blows up. 


3.2 The S^\{Q, Q'}-representation 

The basic conformal structure of the subextremal and extremal Schwarzschild-de Sitter spacetimes 
has already been discussed in laiz] and [47] respectively. Coordinate and Penrose diagrams have 
been also provided in [32] for the subextremal, extremal and hyperextremal cases. In this section 
we present a concise discussion, adapted to our conventions, of the conformal structure of the 
Schwarzschild-de Sitter spacetime in the subextremal, extremal and hyperextremal cases. We 
start our discussion showing that irrespective of the relation of m and A the induced metric at 
the conformal boundary for the Schwarzschild de Sitter spacetime can be identified with the 
standard metric on S^. As discussed in more detail in Section 3.3.1 this construction depends 
on the particular conformal representation being considered. In the subextremal case one cannot 
obtain simultaneously an analytic extension regular near both and re —see [2]. Since we are 
interested only in the asymptotic region, in this section we will consider the region r > r^ For 
the extremal and hyperextremal cases such considerations are not necessary. 

In the following we introduce the null coordinates 




= V + ’^)> 


where r is a tortoise coordinate given by 


r = 


F{r) 


dr. 


(48) 


This integral can be computed explicitly —see mm- The particular form of r depends on the 
relation between A and m. As discussed in [3117] the integration constant can always be chosen 
so that r —0 as r —>■ (X). Defining tanf/ = u,tanl/ = v, with U,V S [—f) f] one gets the line 
element 

1 F'(r) — .,2 TT T/ / Jrr JT/ I 


9SdS = 


sec"^ U sec^ V (dC/ 0 dV + dV 0 dU) — r^cr. 


(49) 
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Figure 4: Penrose diagrams for the extremal Schwarzschild-de Sitter spacetime. Case (a) corre¬ 
sponds to a white hole which evolves towards a de Sitter final state while case (b) is a model of a 
black hole with a future singularity. The continuous black line denotes the conformal boundary; 
the serrated line denotes the location of the singularity; the dashed line shows the location of the 
Killing horizons H at r = 3m . The excluded points Q, Q! and V represent asymptotic regions 
of the spacetime that do not belong to y or the singularity r = 0. 


As discussed in HIE], one can construct Kruskal type coordinates covering the black hole hori¬ 
zon by choosing appropriately the integration constant in equation (48). Analogously, choosing 
a different integration constant, one can construct Kruskal type coordinates covering the cos¬ 
mological horizon. Nevertheless in the subextremal case, as emphasised in [2], it is not possible 
to construct Kruskal type coordinates covering simultaneously both horizons. To construct the 
Penrose diagram for this spacetime, one considers as building blocks the Penrose diagrams for 
the regions 0 < r < rj,, ri, < r < r^. and < r < oo which are then glued together using the cor¬ 
responding Kruskal type coordinates to cross each horizon —see [1132] for a detailed discussion 
on the construction the Penrose diagram and Kruskal type coordinates in the Schwarzschild-de 
Sitter spacetime. Consistent with the above discussion and given that we are only interested 
in the asymptotic region, we restrict our attention, in the subextremal case, to r > Tc. In the 
extremal case one has, however, that r;, = = 3m and one can verify that 


lim 


cos U 


= lim 


cos V 


= C, 


>3m r — 3m r^3m r — 3m 
where C ^ 0 is a constant depending on m and the integration constant chosen in the definition 


of r. Consequently, in the extremal case, the metric (49) is well defined for the whole range of 


the coordinate r: 0 < r < oo —see m- Introducing the coordinates {U, V) defined via 


/ 7r 

tan U = In tan — 
V4 


/tt I/\ 

tan K = In tan —I- 

V4 2 / 


one obtains 


gsds = 2 ■ 


I F{r) 


sec U sec V (d[/ ® dV + dP 0 dU) — 


Recalling that in the subextremal case F{r) < 0 for r > Vc while for the extremal and hyperex¬ 
tremal cases F{r) < 0 for 0 < r < oo, one identihes the conformal factor 


^2 


|A| 


inOI 


cos U cos V. 
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Figure 5: Penrose diagram for the hyperextremal Schwarzschild-de Sitter spacetime. The singu¬ 
larity is of spacelike nature. Dotted lines at 45 ° and 135 ° have been included for visualisation. 
Case (a) corresponds to a white hole which evolves to a final de-Sitter state. Case (b) corresponds 
to a black hole with a future spacelike singularity. 


Therefore, we can identify the conformal metric gsdS = •^‘^dSdS with 

1 - - - - lAIr^ 

OSdS = -X (d[/ 0 dP -I- dP (g) dU) - cost/ cos Per. 

2 |F(r)| 


(50) 


Introducing the coordinates 
one gets 


T=U + V, 


4- = y-t/, 


gsdS = ^ (d^* g) d 4 ' - dr (g) dT) - cos ^ (T -k 4') cos ^ (T - 4') cr. 

The analysis in [5] shows that the conformal factor S tends to zero as r —^ oo. Hence, to identify 
the induced metric at it is sufficient to analyse such limit. Noticing that 


r = 


1 ^ ^ 1 
[V — u) = 




/ tan(7r/4 -k y)A 
\tan(7r/4 -k t/)/ 


and recalling that 


lim r = 0, 

r—^oo 

one concludes that r ^ oo implies 4' = 0 as long as H yk ±^Tr and V yk i^Tr. Using equation 


(431 one can verify that 


|A|r2 


lim I , I 
r-foo |r(r)| 


= 1 . 


Consequently, the induced metric on ^ is given by 


1 T 

h = —dr g) dr — cos^ —cr 
4 2 


which can be written in a more recognisable form introducing ^ | (r -k tt) so that 


h = —d^ g) d^ — sin^ ^(T. 


(51) 
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The metric h is the standard metric on S^. Observe that the excluded points in the discussion of 
this section (C7, V) = correspond to ^ = 0 and ^ = n —the North and South pole of 

S^. The Penrose diagram of the subextremal, extremal and hyperextremal Schwarzschild-de Sitter 
spacetime is given in Figure]^ (a). The conformal boundary of the (subextremal, extremal and 
hyperextremal) Schwarzschild-de Sitter spacetime, defined by the condition S = 0, is spacelike 
consistent with the fact that the Cosmological constant of the spacetime is de Sitter-like —see 
e.g. [HIST]. Moreover, the singularity at r = 0 is of a spacelike nature —see ESI ST]- As pointed 
out in Hina, the Schwarzschild-de Sitter spacetime can be interpreted as the model of a white 
hole singularity towards a final de Sitter state. Alternatively, making use of a reflection 

u I—>■ —u, V I—>■ —V, 

one obtains a model of a black hole with a future singularity —see Figures 

In what follows, we adopt the white hole point of view for the extremal and hyperextremal cases 
so that corresponds to future conformal infinity and we will consider a backward asymptotic 
initial value problem. Consistent with this point of view, for the subextremal case we consider 
asymptotic initial data on and study the development of such data towards the curvature 
singularity located at r = 0 —see Figure 


3.3 The M x §^-representation 


In Section |3.2| we have shown that there exist a conformal representation in which the induced 
metric on the conformal boundary corresponds to the standard metric on S^. A quick inspection 
shows that the metric (511 is conformally flat. In this section we put this observation in a wider 
perspective and show that the induced metric on of a spherically symmetric spacetime with 
spacelike is necessarily conformally flat. In addition, a conformal representation in which 
the induced metric at the conformal boundary corresponds to the standard metric on M x 
is discussed. This conformal representation will be of particular importance in the subsequent 
analysis. 


3.3.1 The conformal boundary of spherically symmetric and asymptotically de Sit¬ 
ter spacetimes 

Following an argument similar to the one given in |43j we have the following construction for 
a spherically symmetric spacetime with spacelike conformal boundary: if a spacetime {Ai,g) is 
spherically symmetric then the metric g can be written in a warped product form 

9 = 1- P^o-, (52) 

where 7 is the 2-metric on the quotient manifold Q = M/SO{3), cr is the standard metric of 
and p : Q — )■ K. If g and g are conformally related, g = Q^g, then the spherical symmetry 
condition for g is translated into the requirement that g can be written in the form 

9 = 1- p'^cr, 

where 7 = ©7 and p = 0p, where 0 does not depend on the coordinates on S^. Near let 
us introduce local coordinates (0,V^) on the quotient manifold Q = M/SO{3) so that 0 = 0 
denotes the locus of Since the conformal boundary is spacelike we have that g(d0, d0) > 0. 
Therefore, the metric induced on by g has the form 

h = —A('0)d'0 (g) dip — p'^(tp)cr, 

where A^tp) is a positive function. Redefining the coordinate ip we can rewrite h as 

h = —p'^{tp){df} (g d^ -I- cr). 
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It can be readily verified —say, by calculating the Cotton tensor of h — that the metric h is 
conformally flat. In Section [3^ it will be shown that, in view of the conformal freedom of the 
setting, a convenient choice is to consider a conformal representation in which the the 3-metric 
on ^ is given by 

h = —dip (g) dip — cr. (53) 


This metric is the standard metric of the cylinder M x with ip G (— 00 , 00 ). It can be verihed 
that this conformal representation is related to the one discussed in Section 3.2 via h = uj'^h, 
where the conformal factor uj and the relation between the coordinates are given by 


V'(C) ='0*-ln|csc^-|-cot^|, w(^) =csc(^). 


(54) 


Equivalently, one has that 


^{ip) = arccos 


g2(v>*-b) — l\ 
g2(y'*-!/') -I- 1 y ’ 


uj{ip) 


2 e^ 


(g 2 b* 




where '0* is a constant of integration. We can directly observe that in this representation ^ = 0 
and ^ = TT correspond to ip = —00 and ip = 00 , respectively. 


3.3.2 The extrinsic curvature of the conformal boundary in the MxS^ representation 


A particularly simple conformal representation for the Schwarzschild-de Sitter spacetime can be 
obtained using the discussion of Section [3.3.1[ Accordingly, take the metric of the Schwarzschild- 
de Sitter spacetime as written in equation (42) with F{r) as given by the relation (43) and 
consider the conformal factor 5 = 1/r. Introducing the coordinates q = Ijr and 0 = 
the conformal metric 


9 = ^'^OeSdS 


is given by 


9 = 


|A| 


(^g'^ - 2mg^ - ^ |A|) dC (g) dC - - 2mg^ - ^|A|) 


-1 


d^i (g) dp — cr. 


The induced metric on the hypersurface described by the condition .n = 0 is given by 

h = —d0 (g d0 — cr. 


It can be verified that g satisfies a conformal gauge for which the conformal boundary has 
vanishing extrinsic curvature. To see this, consider a g-orthonormal coframe {tu“} with 

- 2mp^ - i|A|^ dC, 

and a cr-orthonormal coframe. Denote by {6^} the corresponding dual frame. Using 

this frame we can directly compute the Friedrich scalar s = |V^VcS -|- S —see Appendix 
[D| The computation of the Ricci scalar yields 

R = -I2mg. (55) 


A direct calculation using 


1 


V-detg 

shows that VaV“S = 6mp^ — 2g. Consequently, the scalar s vanishes at the hypersurface defined 
by 2 = p = 0. Contrasting this result with the solution to the conformal constraints given in 


equations (34a)-(34b) we conclude that in this representation the hypersurface described by = 0 


has vanishing extrinsic curvature as claimed. 
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Remark 4. Notice that, in this representation the curvature singularity, located r = 0, corre¬ 
sponds to g = oo. Consequently, is at an infinite distance from the conformal boundary. 

Observe that, the components of the Weyl tensor with respect to the orthonormal frame {eg,} 
as described above are given by 

^1212 = —‘img, Ciaia = mg, Cioio = C2323 = wPi ^2020 = C'soao = 2mp. 

This informat ion w ill be required in the discussion of the initial data for the rescaled Weyl tensor 
—see Section [ 3 . 4. 2 1 Using now that dated = ^~^Cabcd with S = ^ and exploiting the fact that 
the computations have been carried out in an orthonormal frame so that C^'ted = v’^^Cftcdi we 
get 

di2i2 = —2m, di3i3 = m, dioio = —m, ^2323 = w, (^2020 = —m, (i3030 = 2m. 

Finally, considering dij = diojo we have 

dll = —m, d22 = —m, ^33 = 2m. (56) 


3.4 Identifying asymptotic regular data 


As discussed in Section [3.1[ there is a conformal representation in which the induced metric on the 
conformal boundary of the Schwarzschild-de Sitter is the standard metric ft on S^. Nevertheless, 
the asymptotic points Q and Q', as depicted in the Penrose diagram of Figure]^ are associated 
to the behaviour of those timelike geodesics which never cross the horizon —see Appendix 
Despite that, from the point of view of the intrinsic geometry of these asymptotic regions 
—corresponding to the North and South poles of — are regular, from a spacetime point of 
view they are not. This issue will be further discussed Section [3A^ where it will be shown that 
the initial data for the electric part of rescaled Weyl tensor is singular at Q and Q'. Fortunately, 
as exposed in Section [2.4.2| one can exploit the inherent conformal freedom of the setting to select 
any representative of the conformal class [ft] to construct a solution to the conformal constraint 
equations. Taking into account the previous remarks it will be convenient to choose the conformal 
representation discussed in Section 3.3 ft = w^ft with w and ft given in equations (53) and (541, 


in which the points Q and Q' are at infinity respect to the metric ft. 

3.4.1 A frame for the induced metric at 

Consistent with the discussion of the last section, on one considers an adapted frame {I, m, fh} 


such that the metric (53) can be written in the form 

h = —{I + cr) 

where ^ 

I = dt/j, a =-{rn ^ fh + fn ^ m). 

In terms of abstract index notation we have 

hj^j = — Ij^lj — ‘2m(^^mjy 

The frame {I, m,fh} satisfies the pairings 

IjP = —1, mjfhP = —1, Ijm^ = Ijfh^ = mjm^ = fhjfhP = 0. 


(57) 


(58) 
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3.4.2 Initial data for the rescaled Weyl tensor 


The procedure for the construction of a solution to the conformal constraints at the conformal 
boundary requires, in particular, a solution to the divergence equation (35) for the electric part 
of the rescaled Weyl tensor. The requirement of spherical symmetry of the spacetime can be 
succinctly incorporated using the results in [44]. If the unphysical spacetime {M,g) possesses a 
Killing vector X then the initial data encoded in the symmetric tracefree tensor dij must satisfy 
the condition 

£xd^j = 0, (59) 


where £x denotes the Lie derivative in the direction of X on the initial hypersurface. The only 
symmetric tracefree tensor dij compatible with the above requirement is given by 


dij — T hij'j. 


(60) 


where c = dij^P. 

TT-tensors on K^. The general form of symmetric, tracefree and divergence-free tensors (i.e. 
TT-tensors) in a conformally flat setting are well-known —see e.g. |3Jin]. For convenience of 
the reader, in this short paragraph, we adapt the conventions and discussion given in the latter 
references to the present setting. The general the solutions to the equation 

b^dij = 0, (61) 

where h = —5 is the flat metric has been given in |^. One can introduce Cartesian coordinates 
(x“) with the origin of located at a fiduciary position O. Additionally, we introduce polar 
coordinates defined via p = The flat metric in these coordinates reads 

h = —dp (Sj dp — p'^cr. (62) 


Using this notation and taking into account the requirement of spherical symmetry encoded in 
equation (59) the flat space counterpart of the required solution is 

d = —^ (Sdp (g) d/9 -I- j , 

where A* is a constant. In order to obtain an analogous solution in conformally related 3-manifolds 
one can exploit the conformal properties of equation (611 using the following: 


ij be a tracefree symmetric solution to D'^dij = 0 where D is the Levi-Civita 


Lemma 6. Let d. 

connection ofh. Let h = uj^h, then dij = uj 
where D is the Levi-Civita connection ofh. 


^dij is a symmetric tracefree solution to D’‘dij = 0 


This lemma can be found in j^. Here we have adapted the statement to agree with the 
conventions of this article. 

TT-tensors on and K x S^. One can exploit Lemma to derive spherically symmetric 
solutions of the divergence equation (61) in conformally flat 3-manifolds. In particular, the 
metrics fi and h as given in equations (511 and (62) are related via 

h = uj'^h. 


where 

P(0 = cot(^/2), a;(^) = 2sin2(^/2), (63) 

The coordinate transformation p{f) corresponds to the stereographic projection in which the 
origin O of is mapped to the South pole on S^. Alternatively, one can also derive 

p(^) = tan(^/2), uj{f) = 2cos^{£,/2), (64) 


27 







corresponding to the stereographic projection in which the origin of 1 
pole of S^. Using Lemma with equations (631 or (64) one obtains 


d = 




2^1-uj^O 


(3d^ 0 + ft). 


is mapped to the North 


(65) 


Observe thatdU is singular when uj{^) = 1 which corresponds to ^ = 0 and ^ = tt according 
to equations (63) and (64), respectively. Therefore, in this conformal representation the electric 


part of the rescaled Weyl tensor is singular at the North and South poles of S^. Proceeding in a 
analogous way as in the previous paragraphs one can observe that the metrics h and h given in 


equations (53) and (62) are related via 


h = Lo^h 


where 


P(V') = 


- 




- 


A straightforward computation using Lemma renders 

d = Ai, (3d'0 (g) d-^ + ft). 


( 66 ) 


Moreover, since it follows that verifying that dij satishes the condition (59) 


reduces to the computation of Wi = £xh and showing that the components of uji along any leg 
of the frame vanishes —that is 


/'u;,; = 0, 


rrdiOi = 0 , 


rrdujj = 0. 


The latter can easily be done using the Killing equation £xhij = 2D(iX^) = 0 along with 
equations (57) and (58). Finally, comparing expression (66) with equation (56) we can recognise 


that A* = m. Observe that this identification is irrespective of the extrinsic curvature of 


3.5 Asymptotic initial data for the Schwarzschild-de Sitter spacetime 

In the last section it was shown that the K x S^-conformal representation leads to regular asymp¬ 
totic data for the rescaled Weyl tensor. In this section we complete the discussion the asymptotic 
initial data for the Schwarzschild-de Sitter spacetime in this conformal representation. To do so, 
we make use of the procedure to solve the conformal constraints at the conformal boundary as 
discussed in Section |2.3.4| and the specihc properties of the Schwarzschild-de Sitter spacetime. 


3.5.1 Initial data for the Schouten tensor 

Computing the Schouten tensor Sch[h] of ft we get that 

Sch[h] = ® + 2 ^’ 


Equivalently, in abstract index notation one writes 


1 


hj — 


Thus, recalling the solution to the conformal constraints given in equation (34b) we get, 

Lij = —lilj — -(1 — K )hij. 
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3.5.2 Initial data for the connection coefficients 

In order to compute the connection coefficients associated with the coframe {uji} recall that 
= d'f/j and are cr-orthonormal. Equivalently, one has that {ej} = {9,/,, 61 , 62 } with 


ei 


1 

71 


(m + m), 


where (T = m®m + m® m, so that 


62 



m). 


h = (g) 


The connection coefficients can be obtained using the first structure equation (116a) given in 
Appendix |C.l[ Proceeding in this manner, by a straightforward computation, one can show that 
the only non-zero connection coefficient is 72 ^ 1 - In terms of the Ricci-rotation coefficients, the 
latter corresponds to 2-\/2Re(a*) where a* = —^rh’^Snia in the standard NP notation —see m- 
Therefore, the only no-trivial initial data for the connection coefficients is 


72^1 = ^2(0* -I- a*). 


Remark 4. The frame over the cylinder K x introduced in this section is not a global one. 
Nevertheless, it is possible to construct an atlas covering M x such that one each of the charts 
one has a well defined frame of the required form. 

3.5.3 Spinorial initial data 

In this section we discuss the spinorial counterpart of the asymptotic initial data computed in 
the previous sections. 


3.5.4 Spin connection coefficients 

The spinorial counterpart of the asymptotic initial data constructed in the previous sections is 
readily obtained by suitable contraction with the spatial Infeld-van der Waerden symbols —see 
Appendix C.3 Following the discussion of Section 3.5.2 let uj^ = dip and let denote 

an (T-orthonormal coframe. Using equations (123b) of Appendix C.3 we have that the spinorial 
coframe is given by 


, ,AB ^ AB, 

(jJ — Ui 


= {^yAB ^^AB-^^1 ^-^yAB _ 


^AB 


)a ;2 - 


(67) 


Alternatively, one has that the spinorial frame is given by 


eAB = XAB^a 


+ V^UABey'^m^ 


V2zABez mf 


where e^,^. By 


y~ denote the only non-vanishing frame coefficients. Equation (67) allow us 
to compute the reduced connection coefficients ja^cd using the first Cartan structure equation 
(122a) in Appendix C.3 Alternatively, one can use the results of Section 3.5.2 and the spatial 
Infeld-van der Waerden symbols to compute 


CD — j ^ i^CD ^ k 
lAB EF = li kCTAB Cf jCTEF , 


where 

li^k = 


with 

72^2 =-'/2(a* + d*), 72^1 = ^ 2 ( 0 *-h d*). 
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Using the identities ( 123al-( fl23b[ ) in Appendix C.3 one obtains 

lAB^^ EF = 2-\/2(a* + ai,){yAB — ZAB){yEFZ^^ — y^^ Zef)- 
Thus, the reduced connection coefficients are given by 


E) _ ^ CD \ D t \ 

lAB F = 2 ^AB CD = (a* + F[yAB — Zab)- 


( 68 ) 


By computing the spinor version of the connection form f = Iab^f^'^^ using equations 


(68) and (67) one can readily verify that the first structure equation is satisfied. Additionally, 


using the reality conditions. 


XAB^ = —XAB, 


yAB^ = zab, 


ZAB^ = yAB 


we can verify that 'Jabcd is an imaginary spinor —as is to be expected from the space spinor 
formalism. The field ^abcd represents the initial data for the field ^abcd —the imaginary 
part of the reduced connection coefficient Tabcd- The real part of Tabcd corresponds to the 


Weingarten spinor xabcd which, in accordance with equation (34a), is given initially by 

XABCD = KhABCD- 


Rewriting the reduced connection coefficients (68) in terms of the basic valence-4 spinors intro¬ 


duced in Section 4.1 we get for £,abcd — Xabcd the explicit expression 


£,ABCD = -{<- 


a*)(e^ABCD + abcd) 
1 


^*)^Ac{yBD + zbd) + -I- ai,)eBDiyAc + Zac)- 


2V2 


2V2 


3.5.5 Spinorial counterpart of the Schouten tensor 

The spi norial counterpart of the Schouten tensor lij can be directly read from the expressions in 
Section 


3.5.1 


Observe that the elementary spinor x^^ corresponds to the components of k with 


respect to the coframe (67) since 


Ui^^XAB = —X^^XabI^^ = 


dtp = 1 . 


Replacing hij by its space spinor counterpart Habcd we obtain 

hi IaBCD = —XabXcD — Ti^ABCD- 


Equivalently, recalling that the space spinor counterpart of the tracefree part of a tensor l{ij} = 
hj ~ corresponds to the totally symmetric spinor I(abcd) h follows then from 


that 


hj — ^{ij} , 


IaBCD — l-(ABCD) + ^IhABCD- 


Thus, using that for the metric (53) one has r = —2 and that I = = jx, it follows that 


I — — \ and I(abcd) — —X(abXcd) 


= —‘^h^BCD- Therefore, we get 


1 


Iabcd = —2e abcd — -xhABCD- 

6 


(69) 


Finally, recalling the expressions for the components of the spacetime Schouten tensor given in 


(34b) we conclude 


LaBCD = ABCD — g(l — iK^)hABCD- 
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3.5.6 Initial data for the rescaled Weyl spinor 

Following the approach employed in last section, the spinorial counterpart of (661 is given by 


dABCD — Ai,{3lABlcn + hABCo)- 


However, the trace-freeness condition simplifies the last expression since (Ti = 0 implies that 
dij = d[ijy. Therefore dABCD — d(^ABCD) = ^^*1{abIcd)- As the elementary spinor xab can 
be associated to the components of I respect to the co-frame (671 one gets that 


dABCD = ?>A-^X(^abXcD)- 


This last expression can be equivalently written in terms of the basic valence-4 space spinors of 
Section 14.11 as 

(t'ABCD = drne^ ABC D- 

where, in the absence of a magnetic part, we have identified 4>abcd initially with dABCD- 
Observe that have set = m consistent with the discussion of Section |3.4.2[ 


4 The solution to the asymptotic initial value problem for 
the Schwarzschild-de Sitter spacetime and perturbations 

As already discussed in the introductory section, recasting explicitly the Schwarzschild-de Sitter 


spacetime as a solution to the system of conformal evolution equations (29a)-(29i I requires solving. 


in an explicit manner, the conformal geodesic equations. This, as discussed in Appendix |A.2[ 
is not possible in general. Instead, an alternative approach is to study directly the conformal 


evolution equations (29a)-(29i) making explicit the spherical symmetry of the solution and the 
asymptotic initial data corresponding to the Schwarzschild-de Sitter spacetime. This approach 
does not only extract the required information about the reference solution —in the conformal 
Gaussian gauge— but, in addition, is a model for the general structure of the conformal evolution 
equations. The relevant analysis is discussed in Sections |4.1| and [4^ As a complementary analysis, 
we study the the formation of singularities in the evolution equations. In order to have a more 
compact discussion leading to the Main Result, the analysis of the formation of singularities is 
presented in Appendix Finally, in Section |4.3[ we use the theory of symmetric hyperbolic 
systems contained in m to obtain a existence and stability result for the development of small 
perturbations to the asymptotic initial data of the Schwarzschild-de Sitter spacetime. 


4.1 The spherically symmetric evolution equations 

Hitherto, the discussion of the extended conformal Einstein field equations and the conformal 
constraint equations has been completely general. Since we are interested in analysing the 
Schwarzschild-de Sitter spacetime as a solution to the conformal field equations one has to incorpo¬ 
rate specific properties of this spacetime. The most important assumption for our analysis is that 
of the spherical symmetry of the spacetime. Under this assumption, a generalisation of Birkhoff’s 
theorem for vacuum spacetimes with de Sitter-like Cosmological constant shows that the space- 
time must be locally isometric to either the Nariai or the Schwarzschild-de Sitter solutions —see 
m- As the Nariai solution is known to not admit a smooth conformal boundary ElllS], then 
the formulation of an asymptotic initial value problem readily selects the Schwarzschild-de Sitter 
spacetime. 


To incorporate the assumption of spherical symmetry into the conformal field equations en¬ 
coded in the spinorial zero-quantities (19a)-(19d) one has to reexpress the requirement of spherical 
symmetry in terms of the space spinor formalism. In order to ease the presentation we simply 
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introduce a consistent Ansatz for spherical symmetry —a similar approach has been taken in [48j . 
More precisely, we set 


4>ABCD = (t>2 ABCD, 

Qab = xab, 

QaBCD = 02'^ ABCB + 00^^ hABCD, 


^ABCD = Cl ABCD + C 2 ABCD + Cs ABCD + hABCD 

+ ^^{xbb£AC + XaC^Bd) + ^^{yBB^AC + UAC^Bb) 


+ -^izBD£AC + ZaC^Bd), 


XABCD = X2 ^ABCB + gX/i hABCD, 

C AB X AB, ^AB XAB, ^AB VAB, ^AB ZAB, 

Iab = fx xab, 
dAB = dx Xab- 


(70a) 

(70b) 

(70c) 


(70d) 

(70e) 

(70f) 

(70g) 

(70h) 


The elementary spinors xab, Uab, zab, ^abcd hABCD used in the above Ansatz are 
defined in Appendix |C.2[ For further details on the construction of a general spherically symmetric 
Ansatz see inilMl- Alternatively, one can follow a procedure similar to that of Section [3.5.4| — by 
writing a consistent spherically symmetric Ansatz for the orthonormal frame one can identify the 
non-vanishing components of the required tensors. The transition to the spinorial version of such 
Ansatz can be obtained by contracting appropriately with the Infeld-van der Waerden symbols 


taking into account equations (123a)-(123b), (119a)-(119d) and (120a)-( 120c). 


The Ansatz for spherical symmetry encoded in equations (70a)-(70h) combined with the evolu¬ 


tion equations (29a)-(29i) leads, after suitable contraction with the elementary spinors introduced 
in Section |4.1[ to a set of evolution equations for the fields 


(j) 2 , Qx'^, 02 '^, Ca, L, Cz, e°, el, e+,e„, fx- 


This lengthy computation has been carried out using the suite xAct for tensor and spinorial ma¬ 
nipulations in Mathematica —see m- At the end of the day one obtains the following evolution 
equations: 


dre°x = 5X26° - 1x^6° - fx, (71a) 

drel = 5X26° - Ixhcl, (71b) 

= “1X26+ - 5X/ie+, (71c) 

drCf =-^X 2 ef - Ixhef, (71d) 

9 rfx = \x2fx - \xhfx + 0^, (71e) 

drX2 = \xl - lX2Xh - 0f - 0</'2, (71f) 

5rX/« = -5X2 - IXft - 0f> (71g) 

= 1^X26 - 5X/16 - 5X2^y, (71h) 

dr^l = 1^X2?1 - 5X/16 - 5X2^^, (71i) 

dr£.x = - 1 X 2 fx - ^X 2 ^x - ^XhCx, (71j) 

driy = - 1 X 2 ^ 3 + ^X 2 iy - IXh^y, (71k) 

driz = -5X2^1 + ^X 2 iz - 5X/16, (711) 
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= |X20x - \XhQl + 

= gX20| - \XhQt - |x20f + 0</>2, 
dr^l = -gX20| - |x/i0f, 
dr(t>2 = -^X24>2 - Xh(t>2- 


(71m) 

(71n) 

(71o) 

( 71 p) 


The results of the analysis of Sections [3.5.4[|3.53| and [T5?6| provide the asymptotic initial data 
for the above spherically symmetric evolution equations. The resulting expressions are collected 
in the following lemma: 

Lemma 7. There exists a conformal gauge in which asymptotic initial data for the Schwarzschild- 


de Sitter spacetime can he expressed, in terms of the fields defined by the Ansatz (|70a|)-(|70h|), as 
02 =6m, 


ej = 0 , 


Cl = - (a 

e° =0, 

fx=0. 


a*), 

(a* + d*), 


Cs — — (a* + a*), 


X2 =0, 


e' =1, 


02 =- 2 , 

1 


1 


Xh =3/c, 
e'*' =1 


(a* + a*), 


Xx =0, 


4.2 The Schwarzschild-de Sitter spacetime in the conformal Gaussian 
gauge 

In this section we analyse in some detail the spherically symmetric evolution equations derived in 
the previous section. In particular, we show that there is a subsystem of equations that decouples 
from the rest —which we call the core system — and controls the essential dynamics of the system 


(71al-(71p). 


As the Schwarzschild-de Sitter spacetime possess a curvature singularity at r = 0, one expects, 
in general, the conformal evolution equations to develop singularities. Moreover, since the two 
essential parameters appearing in the initial data given in Lemma are m and k —the function 
a* only encodes the connection on — one expects, in general, that the congruence of conformal 
geodesics reaches the curvature singularity at t = (m, k). Nevertheless, numerical evaluations 

suggest that for k = 0 the core system does not develop any singularity —observe that this is 
consistent with the remark made in the discussion of Section [3.3.2[ Furthermore, an estimation 
for the time of existence r@ of the solution to the conformal evolution equations ( |71a[ )- 
with initial data in the case k = 0 is given. A discussion of the mechanism for the formation of 
singularities in the core system (k ^ 0) and the role of the parameter k is given in Appendix 

4.2.1 The core system 

Inspection of the system (7Ia)-(71p) reveals that there is a subsystem of equations that decouple 
from the rest. In the sequel we will refer to these equations as the core system. Defining the helds 

1 

^2'-’ 




x= ■^yipi 2 + xh], 


the system (71p)-(71al can be shown to imply the equations 

0 = -3x0, 

X = -X^ + L- ^00, 


(72) 

(73a) 

(73b) 
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(73c) 


L = -xL- iec;!), 

where the overdot denotes differentiation with respect to r and 


Q{r) = Mr{l + \nr 


The initial data for this system is given by 

(/)(0) = 2m, x(0) = K, 


0 = 


(1 + kt). 




(74) 


As it will be seen in the remainder of this article, equations (73a|-(73cl with initial data (74) 


govern the dynamics of the complete system (71a)-(71p|. The evolution of the remaining fields 


can be understood once the core system has been investigated. 


4.2.2 Analysis of the Core System 

This section will be concerned with an analysis of the initial value problem for the core system 


(73al-( 73c) with initial data given by (74). As it will be seen in the following, the essential feature 


driving the dynamics of the core system (73a)-(73c) is the fact that the function y satisfies a 


Riccati equation coupled to two further fields. One also has the following: 


Observation 1. The core equation (73a) can be formally integrated to yield 


(/>(t) = 2m exp ^—3 J x(s)ds 


(75) 


Hence, (/)(r) > 0 if m ^ 0. 


In the remaining of this section, we analyse the behaviour of the core system in the case where 
the extrinsic curvature of vanishes. 


As discussed in Section 


2.4 


in the case k = 0 the conformal factor reduces to 0 ( 1 ") = ^J\\\/iT 
—thus, one has only one root corresponding to the initial hypersurface J^. To simplify the 
notation recall that 0* = so that 0 (t) = 0*r. Accordingly, the core system (73a)-(73c) 

can be rewritten as 


(j) = -3x((), 

X = -X^ + L- ^0*r(?i, 
L = -xL- 


(76a) 

(76b) 

(76c) 


Moreover, the initial data reduces to 


x(0) = 0, L(0) = 


1 


(/)(0) = 2m. 


Observation 2. A direct inspection shows that equations (76a)-(76c) imply that 

x(t) = tL{t). 


This relation can be easily verified by direct substitution into equations (76b I and (76c). Observe 


that L{t) = x(''‘)/'r is well defined at where r = 0 and x(0) = 0 since the initial conditions 
ensure that 

lim Xll = y 

T->o r 2 
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Taking into account the above observation the core system reduces to 


L = -tL^ - 
(j) = —“irLcj) 


with initial data 



(j){Q) = 2m. 


Observation 3. One can integrate (77b) to 


0 (t) = 2 m exp 



(77a) 

(77b) 

(78) 


(79) 


and conclude that ())(t) > 0 for r > 0. 

To prove the boundedness of the solutions to the core system we begin by proving some basic 
estimates: 


Lemma 8 . If k = 0, then the solution of (73a)-(73c) with initial data (74) satisfies the bound 
L{t) > (fir) ^ 0 . 


Proof. Using equations (77a) and (77b) we obtain the expression 

1 ■ 


(j)L- L(j) = 2TL^(j) - 


for 


T > 0. 


(80) 


Since (j){T) > 0 we can consider the derivative of L/f). Notice that 


dr \(j) 


This observation and inequality (80) gives 


d fL 


^(-1 >- x 0 * 


dr \(j) 


for 


r > 0. 


Integrating the last differential inequality from t = 0 to r > 0 taking into account the initial 
conditions leads to 


L{t) > (j){T) 


f — - ^0*T 

V4m 2 


for 


T > 0. 


□ 


Observe that the last estimate ensures that L{t) is non-negative for r G [O,8m/0*]. It turns 
out that finding an upper bound for L(t) is relatively simple: 


Lemma 9. If k = 0 then, for the solution of (73a)-(73c) with initial data (74), one has that 

2 


L{t) < 


+ 4 


for 


T > 0. 


Proof. Assume r > 0. Using that ^(t) > 0 and equation (77a) one obtains the differential 
inequality 

L{t) < 


35 
















Using that L{t) > 0 for r > 0 one gets 


L{t) 


< —T. 


L^{r) 

The last expression can be integrated giving an upper bound for L{t): 

2 


H't) < 


+ 4’ 


□ 


A simple bound on a finite interval can be found for the field ^(t) as follows: 


Lemma 10. If K = 0 then, for the solution of ( |73a[ )-( [73c| with initial data ( |74[ ) and for 0 < r < 
1/(2 0 *to), the field (j){T) is bounded by above. 

Proof. Assume r > 0. From the estimate of Proposition one has that 

L > 


Therefore 


-3rL(j) < 2 ®*'^ 


2^2 


Using equation (77b) one obtains the differential inequality 

Since ^(t) > 0 the last expression can be integrated to yield, 

2 m 


</(t) < 


1 - 0 * 


Therefore, for 0 < t < 1/ ^ 0*m, the field, (/(r) is bounded by above. Consequently, one can 
take 0 < r < 1 /{2 \/q^^). □ 

The results of Lemmas [HI and [T0| can be summarised in the following: 


Lemma 11. The solution to the core system (73al-(73c) with initial data (74), in the case k = 0, 
is bounded for 0 < r < r,, where 


T, = mm 


r 8m 
0* 2 ^ 0*m 


ft im ^ 


(81) 


Remark 5. A plot of the numerical evaluation of the solutions to the core system (73a)-(73c) 
with initial data (74) in the case k = 0 is shown in Figure]^ 


4.2.3 Behaviour of the remaining fields in the conformal evolution equations 

In this section we complete the analysis of the conformal evolution equations. In particular, we 
show that the dynamics of the whole evolution equations is driven by the core system. To this 
end, we introduce the fields 


X=^{X2-Xh), L=^(0|-0f). 
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Figure 6 : Numerical solution of the core system in the k = 0 case with |A| = 3 and m = l/S-s/S- 
The solid line corresponds to the dashed line to x the dotted line to L. Observe that in 
contrast to the k > 1 and k < — 1 cases, numeric evaluations suggest that in the case k = 0 the 
fields of the core system are bounded for all times —see Figures and IT of Appendix [B| 


The evolution equations for these variables are 

X = -L - 
L = xL + 00, 


(82a) 

(82b) 


with initial data 


1 


L(0) = -d(l + «2). 


x(0) = -K, 

Notice that despite these equations resemble those of the core system, the field 0 is not determined 
by the equations (82a)-(82b) —thus, we call this subsystem the supplementary system. Once the 
core system has been solved, 0 can be regarded as a source term for the system (82a)-(82b|. If 
X and L are known then one can write the remaining unknowns in quadratures. More precisely, 
defining 


Ca = 


1 




^'i/3 0 *^3 


= 6 + 0^1: 




one finds that the equations for these fields can be formally solved to give 

C 3 ('^) = (“ xis)dsj , ^"3 = ^- 3 ( 0 ) exp (^- x(s)ds^ , 

C^i(t) = ?+l( 0 ) exp x{s)dsj , Ci = Ci( 0 ) exp (^- x(s)ds^ . 

The role of the the subsystem formed by 0j, fx and e^ is analysed in the following result. 

Lemma 12. Given asymptotic initial data for the Schwarzschild-de Sitter spacetime, if = 0 
on then 

fxir) = ex°{T) = 0 x ^( t ) = 0 . 
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Proof. This result follows directly from equations (71b),(71e), (71m) and the initial data given 
in LemmaTo see this, first recall that 

doo = x^®eAB*ei( 0 ) = ej;°doQ + 

Assuming then that 63(/c) = 0 one has that 63(8) = 0 and therefore 

dj; = \/2ea;°0. 


Observing that equations (71b),(71e), (71m) form an homogeneous system of equations for the 
fields with vanishing initial data then, using a standard existence and uniqueness 

argument for ordinary differential equations, it follows that the unique solution to this subsystem 
is the trivial solution, namely 

/x(r) = ex°iT) = Qx^{t) = 0 . 


□ 


Using the result of Lemma 12 one can formally integrate equation to yield 

f.x{T) = ^xi.0) exp (^-J x(s)ds^ . 

The frame coefficients can also be found by quadratures 

e°('r) = e°( 0 ) exp x(s)dsj , e+ (r) = e+ ( 0 ) exp (^- x(s)ds^ , 

4 (t) = e+ ( 0 ) exp (^- j x(s)ds^ • 


Since we can write 


X2 =2(x + x), 
1 


~r, (^ 1/3 + 


Xh =2x - X, 

6 =^( 4+1 +Cl), 


0 f = 2 (L-L), 

a =2(a"i-a-i), 


0 f = - Z - 2 L, 
a =2(^+3 - c)- 


then, it only remains to study the behaviour of x and L to completely characterise the evolution 
equations (71a)-(71p). 

Remark 6. In the analysis of the core system of Appendix we identify the mechanism for the 
formation of singularities at finite time in the case k 7 ^ 0. Since (j) acts as a source term for the 
supplementary system (82a)-(82b) one expects the solution to this system to be singular at finite 
time if the solutions to the core system develop a singularity. Clearly, the behaviour of the core 
system is independent from the behaviour of the supplementary system. Consequently, the fact 
that (j) diverges at finite time or not is irrespective of the behaviour of Z and x- 


4.2.4 Deviation equation for the congruence 


As discussed in Section 2.3.2 the evolution equations (29a)-(29h) are derived under the assump¬ 
tion of the existence of a non-intersecting congruence of conformal geodesics. In this section we 
analyse the solutions to the deviation equations. 


As a consequence of Lemma 12 we have Jab = 0. Following the spirit of the space spinor 
formalism, the deviation spinor zab can be written in terms of elementary valence 2 spinors as 


Z{AB) — ZxXaB + ZyUAB + ZzZaB- 
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Substituting expression (70e) into equation (31b I and using the identities given in equation (120d) 
one obtains ^ ^ 

~ bj ~ b, d^-Zy = ^Xh^y 


6 ' 


One can formally integrate these equations to obtain 

Zxir) = z^{T) = z^i,, Zy{T) = Zyi, eyzp J x(s)ds^ . 

In the last equation, Zy* and Zz* denote the initial value of Zx{t), %(t) and Zz(t) respectively. 
It follows that the deviation vector is non-zero and regular as long as the initial data z^,*, Zy* and 
Zz* are non-vanishing and x('^) is regular. Accordingly, the congruence of conformal geodesics 
will be non-intersecting. 

4.2.5 Analysis of the supplementary system 

As in the case of the core system, the supplementary system is simpler in the gauge in which 
K = 0. In such case, direct inspection shows that equations (82a)-(82b) imply 

X = -rL. 


This can be verified by direct substitution into equations (82a) and (82bI. Notice that L{t) is 
well defined at where t — 0 and x(0) = b since the initial conditions ensure that 


ita Ihl 

T-J-O T 


Taking into account this observation, the system (82a)-(82bl reduces to the equation 

t = -tL^ + 0*(/), 

with initial data 

i(o) = -^. 

Using that f) is only determined by the core system, together with the analysis of Section |4.2.2| 
one obtains the following result: 


(83) 


(84) 


y0* -I- 2arctan Q©* 


Lemma 13. The solution to equation (83) with initial data (84) is bounded for 0 <t <t@ with 

r@ = mm{To,T, }, where To = a/0 * 

Proof. To prove that L{t) is bounded from above we proceed by contradiction. Assume that 
L —)■ 00 for some finite S [0 ,t,], then L —> oo at . Now, equation (83) can be rewritten as 

L -\- rlf = 0 *((). 

Therefore, since t > 0, the last expression implies that (/)—>■ oo at . However, in Section 


4.2.2| we showed that cf is finite for r G [b,T,]. This is a contradiction, and one cannot have 
L —>• oo at S [0,T,]. Consequently L(t) is bounded from above for 0 < r < r,. To show 
that L(t) is bounded from below, for 0 < r < To with Tq as given by relation (85), observe that 
(^(t) > —r for r > 0 since f)(T) > 0. Using this observation, equation (83) implies the differential 
inequality 

i>-T(L^ + ej. 
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Since 0* > 0 one has that + 0*) > 0. Thus, one can rewrite the last inequality as 


(T2 + 0*) 


> -T, 


which can be integrated using the initial data (84) to give 

L(r) > — i/^tan ( + arctan ( —| | . 

V2 \2^/¥JJ 

Since the function tan is bounded if its argument lies in [0,7r/4] one concludes that L{t) is 
bounded from below for 0 < r < Tq. Finally, taking the minimum of r, and Tq one obtains the 
result. □ 

Remark 7. Numerical evaluations of the solutions to the supplementary system show that it 
should be possible to improve Lemmas 11 and 13 and conclude that the solutions do not blow up 
in finite time. These results, however, will not be required to formulate the Main Result of this 
article. 


4.3 Perturbations of the Schwarzschild-de Sitter spacetime 


In the sequel, we consider perturbations of the Schwarzschild-de Sitter spacetime which can be 
covered by a congruence of conformal geodesics so that Lemma can be applied. In particular, 
this means that the functional form of the conformal factor is the same for for both the background 
and the perturbed spacetime. 

The discussion of Section [3 . 4| brings to the foreground the difficulties in setting up an asymp¬ 
totic initial value problem for the Schwarzschild-de Sitter spacetime in a representation in which 
the initial hypersurface contains the asymptotic points Q and Q': on the one hand, the initial 
data for the rescaled Weyl tensor is singular at both Q and Q'; and, on the other hand, the curves 
in a congruence of timelike conformal geodesics become asymptotically null as they approach Q 
and Q' —see Appendix]^ 


Consistent with the above remarks, the analysis of the conformal evolution equations (29a)- 


(29h) has been obtained in a conformal representation in which the metric on ^ is the standard 
one on K X . In this particular conformal representation the asymptotic points Q and Q' are at 
infinity respect to the 3-metric of and the initial data for the Schwarzschild-de Sitter spacetime 
is homogeneous. In this section we analyse nonlinear perturbations of the Schwarzschild-de Sitter 
spacetime by means of suitably posed initial value problems. More precisely, we analyse the 
development of perturbed initial data close to that of the Schwarzschild-de Sitter spacetime 
in the above described conformal representation. Then, using the conformal evolution equations 


(29al-(29h) and the theory of first order symmetry hyperbolic systems contained in m we obtain 
a existence and stability result for a reference solution corresponding to the asymptotic region of 
the Schwarzschild-de Sitter spacetime —see Figure 


4.3.1 Perturbations of asymptotic data for the Schwarzschild-de Sitter spacetime 

In what follows, let S denote a 3-dimensional manifold with S « K x S^. By assumption, there 
exists a diffeomorphism '0 : 5 —?► K x which can used to pull-back a coordinate system x = (a;“) 
on K X to obtain a coordinate system on S —i.e. x ='ll}*x = xo'tjj. Exploiting the fact that 
■0 is a diffeomorphism we can define not only the pull-back ip* : T*(M. x S^) —>■ T*S but also the 
push-forward of its inverse {ip~^)* : T(K x S^) — TS. Using this mapping, we can push-forward 
vector fields Cj on T(IR x S^) and pull-back their covector fields a* on T*S via 

Ci = S* = 0*aL 
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In a slight abuse of notation, the fields Cj and 3* will be simply denoted by Cj and a* 


In the following, we will refer to all the helds discussed previously for the exact Schwarzschild- 
de Sitter spacetime as the background solution and distinguish them with a over the Kernel 
letter —e.g. h will denote the standard metric on K x given in equation (531. Similarly, the 
perturbation to the corresponding held will be identihed with a over the Kernel letter. Notice 
that although the frame {cj} is /i-orthonormal, it is not necessarily orthogonal respect to the 
intrinsic 3-metric h on S. 

Let {ej} denote a h--orthonormal frame over TS and let {u)*} be the associate cobasis. Assume 
that there exist vector helds {e^} such that an ft,-orthonormal frame {ej} is related to an h- 
orthonormal frame {cj} through the relation 

— Ci Gj- 

This last requirement is equivalent to introducing coordinates on S such that 

h = h + h. (86) 


Now, consider a solution 


; \.ij ^ ^ijhi ^ij') 


to the asymptotic conformal constraint equations (33a)-(33i| which is, in some sense to be deter¬ 
mined, close to initial data for the Schwarzschild-de Sitter spacetime so that one can write 


— ^ij\s + Xij\s — Xij\s + Xijlsj 

^ij\s ^ij\s “f dijh\s “f ^ijh\si 


Li\s — Li\s Li\g 
^ij 15 ^ij 15 “f ^ij I S • 


A spinorial version of these data can be obtained using the spatial Infeld-van der Waerden 
symbols. Accordingly, one writes 

VabcdIs = VabcdIs + VABcnls, IJ-abcd\s = fJ-ABcnls, (87a) 

LabcdIs = LabcdIs + LabcdIs, ^abcdIs = ^ABcnls + ^abcdIs, (87b) 
Lab\s = Lab\s, XABcnls = XabcdIs + XABcnls, (87c) 

eAB\s = ^Ab\s + ^Ab\s, fAB\s = fAB\s- (87d) 


Observe that all the objects appearing in expressions (87a)-(87dl are scalars. 


4.3.2 Controlling the size of the perturbation 

In this subsection we introduce the necessary notions and dehnitions to measure the size of the 
perturbation of the initial data. Let A = {{(j)i,Ui), (<(> 2 ,^ 2 )} with (j)i :Ui ^ and ^2 : ^2 —t 
be an Atlas for K x S^. Let Vi C V 2 C U 2 be closed sets such that K x C Vi U V 2 . In 
addition, define the functions 


Vi{x) = 


1 xe(/>i(Vi) 


0 X £ 


(Vi) 


rj2{x) = 


X G (('2(V2) 


0 a; G 


i(V2) 


( 88 ) 


Observe that any point p G S is described in local coordinates by Xp = (^i o'0)(p) with Xp G <(>(77) 
where ip is the diffeomorphism defined in Section 4.3.1 and {(pM) G A. Consequently, any smooth 
function Q : 3 —>■ can be regarded in local coordinates as Q{x) : (p{U) -A . Let Qi{x) 

denote the restriction of Q{x) to one the open sets (piiUi) for i = 1, 2. Then, we define the norm 
of Q as 

I! Q l|5,m = || 'ni{x)Ql{x) ||R3_m -h II f] 2 {x)Q 2 {x) ||r3_„ 
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where 


Q IIb3 


C m 3 p ^ 

y] y] / {da,--da,Qfd^x 

1=0 o:i,...,ai > 


1/2 


Now, we use these notions to define Sobolev norms for any quantity Qic with k: being an 
arbitrary string of frame spinor indices as 


Qk. I!s ,m — E II Qk ||(S,m ■ 


In the last expression m is a positive integer and the sum is carried over all the independent 
components of Qtc which have been denoted by Q^. 

4.3.3 Formulation of the evolution problems 


Consistent with the split (87a)-(87dl for the initial data, we look for solutions to the conformal 
evolution equations (38a)-(38b) of the form 

VABCD = VABCD + f]ABCD, fJ-ABCD — i^ABCD, (89a) 

LaBCD = LaBCD + LabCB, ^ABCD = ^ABCD + ^ABCn, (89b) 

Lab = Lab, Xabcd = Xabcd + Xabcd, (89c) 

Gab = Cab + cab, Jab = Jab- (89d) 


Using the notation introduced in Section 2.5 the initial data (87a)-(87d) will be represented 
as u*. The perturbed initial data will be assumed to be small in the sense that given some e > 0 
one has 


lltS.m— II \A.HCD ||(S,m + II ^ABCn lls 
+ II Cab lls ,m + II /as lls ,m + 


II LabCD lls ,m H" 
iABCD ||s,m< £■ 


Lab 


|»S,m 


Remark 8. Notice that, as a consequence of the conformal representation being considered, the 
above smallness requirement on the perturbed initial data constraints the possible behaviour of 
the perturbation near the asymptotic points Q and Q!. To see this in more detail let 0 denote 
a perturbation of the initial data for some component the rescaled Weyl spinor. For simplicity, 
assume that in some local coordinates ('0,6*,(/)) for M x S^, the perturbed field (j) is independent 
of In such case, ii ^ G L^(M) one has that 

^ (90) 


with f3 > 1/2. Consequently, in the K x S^-conformal representation the perturbations must 
decay at infinity, i.e. as they approach Q and Q'. Under the conformal transformation g = w'^g 
the components of the rescaled Weyl spinor transform as 4>abce> = ^~^4>abcd- This last 
expression is consistent with the frame version of the conformal transformation rule given in 
Lemma 1^ Taking into account the discussion of Section 3.3.1 and equation (90) one concludes 
that for the corresponding perturbation in the S^-conformal representation one has 

<^=o(r'(iniei)-^) 


near the South pole ^ = 0. Consequently, initial data on K x satisfying L^-decay conditions 
near infinity correspond, in general, to data which is singular in other conformal representations. 
In other words, the class of perturbation data that we can consider can be, in principle, singular 
at both the North and South poles in the S^-conformal representation. 
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(a) 


(b) 



Figure 7: Schematic depiction of the development of perturbed initial data for the Schwarzschild- 
de Sitter spacetime and the congruence of conformal geodesics. In (a) the evolution of asymptotic 
initial data is depicted in the conformal representation in which the asymptotic points Q and Q' 
are at a finite distance respect to the metric on Figure (b) shows a schematic depiction of the 
evolution of asymptotic initial data in the conformal representation in which Theorem has been 
formulated. In contrast to the conformal representation leading to Figure (a) , the initial data 
is homogeneous and formally identical for the subextremal, extremal or hyperextremal cases. In 
both diagrams, the dashed line corresponds to the location of an hypothetical Cauchy horizon of 
the development. 


Remark 9. An explicit class of perturbed asymptotic initial data sets can be constructed, keeping 
the initial metric fixed to be standard one on R x S^, using the analysis of [5] as follows: introduce 
Cartesian coordinates (a:“) in R^ with origin located at a fiduciary position Q and define a polar 
coordinate via p = Sapx'^x^. The general solution of the equation 


D^d,j = 0 , 


where £)* is the Levi-Civita connection on R^, can be parametrised as 

liP) ^ ^ ^ i(Q) ^ d(A) 


da, = d'-:>+d>->+d 


^ah ' 


The terms are divergent at Q and have been explicitly derived in [5] . Given 

any smooth function A(x“) on R^ the term can be obtained using the operators 3 and 8 —see 
m for definitions. This term can have, in general, any behaviour near Q —see [3]. However, 
setting A = 0{p'^) with n > 3 the term is regular near Q. Using the frame version of the 
conformal transformation rule of Lemma|^and either equation (|63|) or (64) one can verify that the 


corresponding term in the S^-representat ion is 4^) = C>(p’^+3). Similarly, using the conformal 
transformation formulae, given in Section [3.3.l| relating the and R x S^-representations of the 
initial data, one obtains d!'^ = We observe that regular behaviour of perturbed initial 

data in the R x S^-representation does not necessarily correspond to regular behaviour in the 
S^-representation nor in the R^-representation. 


4.4 The main result 

The main analysis of the background solution in Section [4.2| was performed in a conformal repre¬ 
sentation in which the asymptotic initial data is homogeneous and the extrinsic curvature of 


vanishes —i.e. k = 0. The general evolution equations (38al-(38b) consist of transport equations 


for V coupled with a system of partial differential equations for cf). However, as shown in Section 
4.2 the assumption of spherical symmetry implies that the only independent component of the 


spinorial field 4>abcl> is (/> 2 - Consequently, the system (38a)-(38bl reduces, for the background 


fields u = (v,(j)), to a system of ordinary differential equations. The Piccard-Lindeldf theorem 
can be applied to discuss local existence of the latter system. However, one does not have, a 
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priori^ control on the smallness of the existence time. To obtain statements concerning the ex¬ 
istence time of the perturbed solution, we recall that the discussion of the evolution equations of 


Section 4.2 shows that the components of solution u are regular for r G [0,r@] with t@ as given 
in equation (85), so that the guaranteed existence time is not arbitrarily small. 


The analysis of the core system in Section 4.2 was restricted to the case k = 0, in which 
the conformal boundary has vanishing extrinsic curvature. In this case, we obtained an explicit 
existence time r@ for the solution to the conformal evolution equations. In contrast, the analysis 
given in Appendix shows that in general, for k ^ 0, the core system develops a singularity 
at finite . Since the results given in Section 4.2.4 for the conformal deviation equations hold 


not only for k = 0 , but for any k as long as = 0 , one has that the congruence of conformal 
geodesics is non-intersecting in the k ^ 0 case as well. This shows that, the singularities in the 
core system in the case k 7 ^ 0 are not gauge singularities. The estimation for the existence time 
T@ in the k = 0 case along with the discussion of the reparametrisation of conformal geodesics 
given in Appendix |B.3| can, in principle, be used to obtain an estimation for the existence time 
T 0 in the case k 7 ^ 0 . 

In this section it is shown how one can exploit these observations, together with the theory 
for symmetric hyperbolic systems, to prove the existence of solutions to the general conformal 
evolution equations with the same existence time t@ for small perturbations of asymptotic initial 
data close to that of the Schwarzschild-de Sitter reference solution. By construction, the devel¬ 
opment of this perturbed data will be contained in the domain of influence which corresponds, 
in this case, to the asymptotic region of the spacetime —see Figure 

Taking into account the above remarks and using the theory of symmetric hyperbolic systems 
contained in m one can formulate the following existence and Cauchy stability result: 

Theorem 1 {existence and Cauchy stability for perturbations of asymptotic initial 
data for the Schwarzschild-de Sitter spacetime). Let u* = -|- denote asymptotic 

initial data for the extended conformal Einstein field equations on a 3-dimensional manifold S « 
K X where u* denotes the asymptotic initial data for the Schwarzschild-de Sitter spacetime 
(subextremal, extremal and hyperextremal cases) with k = 0 in which the asymptotic points Q and 
Q! are at infinity. Then, for m > A and r@ as given in equation (85), there exists £ > 0 such 
that: 


(i) for ||u*|| 5 „i < £ , there exist a unique solution ii to the conformal evolution equations 
with a minimal existence interval [ 0 , r@] and 

uGC™-2([0,r@] x5,C^), 

and the associated congruence of conformal geodesics contains no conjugate points in [ 0 ,r@]; 

(ii) given a sequence of perturbed data such that 

II ||s,m-t 0 as 00 , 

then the corresponding solutions have a minimum existence interval [0 ,t@] and it 

holds that 

II lU.m-t 0 as 00 

uniformly in t G [0, t®] as n ^ 00 ; 

(Hi) the solution u = u-|-u is unique in [0,r@] x S and implies a solution (A4r@;9) to 

the Einstein vacuum equations with the same de Sitter-like Cosmological constant as the 
background solution where 

Mr. = (0,T@) X S. 

Moreover, the hypersurface = {0}x5 represents the conformal boundary of the spacetime. 


(40a)-(40b 


44 
















Proof. Points (i) and (ii) are a direct application of the theory contained in |37j where it is used 
that the background solution u is regular on t G The initial data for the Schwarzschild- 

de Sitter spacetime encoded in u* is in a representation in which the points Q and Q' are at 
infinity. Observe that the asymptotic initial data, as derived in Section [3.5[ for the subextremal, 
extremal and hyperextremal cases are formally the same —in particular, notice that the initial 
data for the electric part of the rescaled Weyl tensor contains information about the mass m 
while the conformal factor 0 carries information about A. The arguments in the analysis of 
Section |4.2| are irrespective of the relation between A and m. The key observation in the proof is 
that one can apply the general theory of symmetric hyperbolic systems of m for each open set 
and chart of an atlas for K x S^; then, these local solutions can be patched together to obtain 
the required global solution over [0, r@] x S —it is sufficient to cover K x with finitely many 
patches (two) as discussed in Section 4.3.2 Details of a similar construction in the context of 
characteristic problems can be found in m To prove point (in) first observe that from Lemma 

-(40b I implies a solution u = u -|- u to the 
X 5 if u* = u* -|- u* solves the conformal 


[^the solution to the conformal evolution system (|40a 


extended conformal Einstein field equations on [0, t@ 
constraint equations on the initial hypersurface. This solution implies, using Lemma[^ a solution 
to the Einstein field equations whenever the conformal factor is not vanishing. General results of 
the theory of asymptotics implies then that the initial hypersurface S can be interpreted as the 
conformal boundary of the physical spacetime —see [5TJ|331. □ 


5 Conclusions 

In this article we have studied the Schwarzschild-de Sitter family of spacetimes as a solution to the 
extended conformal Einstein field equations expressed in terms of a conformal Gaussian system. 
Given that, in principle, it is not possible to explicitly express the spacetimes in this gauge, we 
have adopted the alternative strategy of formulating an asymptotic initial value problem for a 
spherically symmetric spacetime with a de Sitter-like Gosmological constant. The generalisation 
of Birkhoff’s theorem to vacuum spacetimes with Cosmological constant then ensures that the 
resulting solutions are necessarily a member of the Schwarzschild-de Sitter spacetime. 

As part of the formulation of an asymptotic initial value problem for the Schwarzschild-de 
Sitter spacetime we needed to specify suitable initial data for the conformal evolution equations. 
The rather simple form that the conformal constraint equations acquire in the framework con¬ 
sidered in this article allows to study in detail the conformal properties of the Schwarzschild-de 
Sitter spacetime at the conformal boundary and, in particular, at the asymptotic points where 
the conformal boundary meets the horizons. The key observation from this analysis is that the 
conformal structure is singular at these points and cannot be regularised in an obvious manner. 
Accordingly, any satisfactory formulation of the asymptotic initial value problem will exclude 
these points. 

An interesting property of the conformal evolution equations under the assumption of spherical 
symmetry is that the system reduces to a set of transport equations along the conformal geodesics 
covering the spacetime. The essential dynamics, and in particular the formation of singularities 
in the solutions to this system, is governed by a core system of three equations —one of them 
a Riccati equation. As discussed in Appendix this core system provides a mechanism for 
the formation of singularities in the exact solution. The analysis of the core system allows 
not only to study the properties on the Schwarzschild-de Sitter spacetime expressed in terms 
of a conformal Gaussian gauge system, but also to understand the effects that the gauge data 
has on the properties of the conformal representation arising as a solution to the conformal 
evolution equations. It is of interest to explore the idea of whether the mechanisms identified in 
the analysis of the core system could be used to analyse the formation of singularities in more 
complicated spacetimes —say, in the developments of perturbations of asymptotic initial data for 
the Schwarzschild-de Sitter spacetime. 


45 










The conformal representation of the Schwarzschild-de Sitter spacetime obtained in this article 
has been used to show that it is possible to construct, say, future asymptotically de Sitter solutions 
to the Einstein vacuum Einstein with a minimum existence time —as measured by the proper 
time of the conformal geodesics used to construct the gauge system— which can be understood 
as perturbations of a member of the Schwarzschild-de Sitter family of spacetimes. As already 
mentioned in the main text, it is an interesting problem to determine the maximal Cauchy 
development to these spacetimes. In order to obtain the maximal Cauchy development of suitably 
small perturbations of asymptotic data for the Schwarzschild-de Sitter one would require the use of 
more refined methods of the theory of hyperbolic partial differential equations as one is, basically, 
confronted with global existence problem for the conformal evolution equations. In this respect, 
we conjecture that the time symmetric conformal representation in which k = 0 together with the 
global stability methods of [3H] should allow us to make inroads into this issue. Closely related to 
the construction of the maximal development of perturbations of asymptotic initial data of the 
Schwarzschild-de Sitter spacetime is the question whether there is a Cauchy horizon associated to 
the boundary of this development. If this is the case, one would like to investigate the properties 
of this horizon. Intuitively, the answer to these issues should depend on the relation between 
the asymptotic points Q and Q' and the conformal structure of the spacetime. In particular, 
one would like to know whether the singularities of the rescaled Weyl tensor at these points 
generically propagate along the boundary of the perturbed solution —notice, that they do not 
for the background solution. If one were able to use the M x S^-representation of the conformal 
boundary of perturbations of asymptotic initial data for the Schwarzschild-de Sitter to construct 
a maximal development and to gain sufficient control on the asymptotic behaviour of the various 
conformal fields, one could then rescale this solution to obtain a representation with a conformal 
boundary of the form \ {Q, Q'}. As discussed in the main text, in this representation some 
fields are singular at Q and Q'. This observation suggests that this construction could shed some 
light regarding the propagation (or lack thereof) of singularities near the asymptotic points Q 
and Q'. 
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A Appendix: The asymptotic points Q and Q! and confor¬ 
mal geodesics in the Schwarzschild-de Sitter spacetime 

A.l Analysis of the asymptotic points Q and Q' 

In Section [372] it was shown that there exist a conformal representation of the Schwarzschild-de Sit¬ 
ter spacetime in which the metric at the conformal boundary is h —i.e. the standard metric on S^. 
In addition, we observed that the North and South pole of correspond to special points in the 
conformal structure that we have labelled as Q and Q'. These asymptotic regions are represented 
in the Penrose diagram for the subextremal, extremal and hyperextremal Schwarzschild-de Sitter 
spacetime as the points where the conformal boundary and the Cosmological horizon. Killing 
horizon a nd si ngularity, respectively, seem to meet —see Figures and and As discussed 
in Section |3.2| these points correspond to {U,V) = (±^,±|) for which the tortoise coordinate r 
is not well defined. In Section |3.4| we showed that in the conformal representation in which the 
initial metric is h the data for the electric part of the rescaled Weyl tensordjj, as given in equation 
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(65), is singular precisely at Q and Q'. Observe that written in spinorial terms the initial data 
for the rescaled Weyl spinor in this conformal representation is given by 


<pABCD = 


6 m 




e ABCD 


which is singular at both Q and Q'. This situation resembles that of the geometry near spacelike 
infinity of the Minkowski spacetime and the construction of the cylinder at infinity given in |22j 
which allows to regularise the data for the rescaled Weyl spinor. However, some experimentation 
reveals that this type of regularisation procedure (in contrast with the analysis of Schwarzschild 
spacetime given in [22] 1 cannot be implemented in the analysis of the Schwarzschild-de Sitter 
spacetime without spoiling the regular behaviour of the conformal factor. Since the hyperbolic 
reduction procedure for the extended conformal Einstein field equations is based on the existence 
of a congruence of conformal geodesics in spacetime, the singular behaviour of the initial data for 
the rescaled Weyl spinor suggest that the congruence of conformal geodesics does not cover the 
region of the spacetime corresponding to Q and Q'. To clarify this point, in the remaining of this 
section we analyse the behaviour of conformal geodesics as they approach the asymptotic points 
Q and Q'. 


A.2 Geodesics in Schwarzschild-de Sitter spacetime 


The method for the hyperbolic reduction for the extended conformal Einstein field equations 
available in the literature requires adapting the gauge to a congruence of conformal geodesics. 
The behaviour of metric geodesics in the Schwarzschild-de Sitter spacetime has been already 
studied [33 [31] and an analysis of conformal geodesics in Schwarzschild-de Sitter and anti-de- 
Sitter spacetimes is carried out in (SD). In static coordinates (t,r,9,(p) the equation for radial 
timelike geodesics, {9 — 9^,, ip = (pi,) with 0* and Pi, constant, are 




dt 7 
dr F(r) 


(91) 


The first equation can be formally integrated as 


r - r* 



1 

- F{s) 


ds 


(92) 


where r is the gsds-proper time and 7 is a constant of motion which can be identified with the 
specific energy of a particle moving along the geodesic. The equation for t can be solved once 
equation (92) has been integrated. As pointed out in [3117], by choosing 7 = 1 one can explicitly 
solve this integral. However in general, for arbitrary 7, the integral is complicated and cannot be 
written in terms of elementary functions. A side observation is that if r ^ and r ^ Vc then the 
curves of constant t correspond to geodesics with 7 = 0 . Finally, its worth noticing that geodesics 
with constant r are characterised by the condition 


72 _ F{r) = 0. 


(93) 


This last type of curves, which will be called critical curves, are analysed in Section |A.4| In 
general, the properties of conformal geodesics differ from their metric counterparts. However, 
in the case of an Einstein spacetime with spacelike conformal boundary any conformal geodesic 
leaving orthogonally is, up to reparametrisation, a metric geodesic — see [33] and Lemma 


A.3 


A special class of conformal geodesics in the Schwarzschild-de Sit¬ 
ter spacetime 


As briefly mentioned in Section A.2 and pointed out in [3117], in general, the integral ( [92| ) cannot 
be written in terms of elementary functions except for the the special case when 7 = 1 where it 
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yields 




(94) 


where C is an integration constant. The last expression is valid irrespective of the relation between 
m and A. One can also use this expression to integrate the second equation in (911 to obtain the 
geodesic parametrised as {r{f),t{f)). The integration of t will not be required for the purposes of 
the analysis of this section. A complete analysis of conformal geodesics in the Schwarzschild-de 
Sitter and anti-de Sitter spacetimes will be given in [30j . By virtue of Lemma one can recast 
the geodesic with 7 = 1 as a conformal geodesic by reparametrising it in terms of the unphysical 
proper time as determined by equations (16 1 and (37). A straightforward computation yields 


7(r) = 


In 


2 + ACT 


(95) 


Equivalently, assuming either n > 0 and r > 0 or ac < 0 and 0 < r < — 2/k one obtains in both 
cases 


r(f) = 


2 exp I \ — 


1 — K exp 

From the last expression one can verify that 


(96) 


lim t(t) = 0, 


lim r(f) = —2 /k, 

r—>-oo 


as expected. In what follows, we will rewrite equation ( |94[ ) in terms of the unphysical proper time 
as 

1 (m|A|C3r3-6(2 + /cT)3)2/3 


™(r) = 


(m|A|)2/3 


From the last expressions one can verify that one has r 
location of the singularity r = 0 is determined by 


Ct{t + 2kt) 

oo as T 


0 and 


(97) 

—21k. The 


2 

“ (m|A|)i/3C-K- 


Recalling that C is an integration constant which depends on the initial data for the congruence, 
since the only freedom left in the conformal factor is encoded in k, one realises that C = C{k). So 
one cannot draw any precise conclusion about the location of the singularity unless one further 
identifies explicitly C{k). In particular, considering constant k as we have done for the analysis of 
the core system and setting C to be proportional to k, say C = for some proportionality 

constant >£■, one obtains 

1 


which is in agreement the with the qualitative behaviour of the core system as shown in Figures 
and[l0l Notice, however, that the arguments of the core system given in Section [4. 2. 2| and 
Appendix |B| do not rely on integrating (92) explicitly as we have done in this section. 


A.4 Critical curves on the Schwarzschild-de Sitter spacetime 

In order to clarify the role of the asymptotic points, in this section we show that there are not 
timelike conformal geodesics reaching Q and Q' orthogonally. More precisely, we show that a 
timelike conformal geodesic becomes asymptotically null as it approaches Q and Q'. This is 
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a) 


b) 


Q 0 ! 



Figure 8 : Curves of constant r and t in the Schwarzschild-de Sitter spacetime, a) Curves with 
constant t and r (red and blue respectively) are plotted on the Penrose diagram of the Subextremal 
Schwarzschild-de Sitter spacetime. Curves of constant t accumulate at the bifurcation spheres B, 
B' while the curves of constant r accumulate at the asymptotic points Q and Q'. b) Curves with 
constant t and r (red and blue respectively) are plotted on the Penrose diagram of the extremal 
Schwarzschild-de Sitter spacetime. In contrast with the subextremal case, curves with constant 
t in starting from some r* < 3m accumulate at the asymptotic points Q and Q' while those 
starting from r* > 3m accumulate at V. The hyperextremal case is qualitatively similar to the 
extremal one and has been omitted. 


in stark tension with the required conditions for constructing a conformal Gaussian system of 
coordinates in the neighbourhood of Q and Q'. 

As shown in the Penrose diagram of Figure, in the subextremal case the curves of constant 
t = ti, accumulate in the bifurcation spheres B and B' while the curves of constant r accumulate 
in the asymptotic points Q and Q'. By contrast, in the extremal case the curves with constant 
t = t-i, approach the asymptotic points Q and Q' —see [32] for an extensive discussion on the 


Penrose diagram for Schwarzschild-de Sitter spacetime. It follows from the geodesic equation (91) 


that the curves of constant r correspond to geodesics whenever the condition (931 is satisfied, 
this equation explicitly reads 

|A|r3-t3(7^-l)r-b6m = 0. (98) 

Observe that for 7 = 1 the last condition reduces to |A|r^ -|- 6 m = 0 which cannot be solved for 
positive r. 

In this section we perform an analysis of the behaviour of the critical curves on the Schwarzschild- 
de Sitter spacetime. Notice that in the hyperextremal case the are no timelike geodesics with 


constant r since for |A| > l/9m^ one has strictly F{r) < 0 so that the condition (93) can never 
be satisfied. 


A.4.1 Critical curves in the extremal Schwarzschild-de Sitter spacetime 

We start the analysis in the simpler case in which |A| = l/9m^ so that Fir) is given as in equation 
(46) and the condition (93) reduces to considering r = 3m and 7 = 0 . Observe that the curves 
with 7 = 0 and r ^ 3m correspond to curves with constant t = t* which, as discussed in previous 
paragraphs, approach asymptotically the points Q and Q'. Notice that for 7 = 0 the expression 


(92) can easily be integrated to yield 

f - T* = 3mln(i7(r)/iJ(r,^)) 


(99) 


where 

H( 1 — + \/r + 6 m 

(•\/3r — ^/r -\- 6m){^/r + Vr + 6 m) 
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Observe that equation ( |^ , as pointed out in [47], implies that the geodesics with 7 = 0 never 
cross the horizon since r —>■ 00 as r —>■ 3 to. For simplicity, let AF*- = + exp(T*/3m) with 

r* ^ 3m so that r = 3mIn |i7(r)/M*|. Reparametrising using equation ( j^ and that |A| = l/9m^ 
renders 

ill -iw(r) 

with W{r) = Using L’Hopital rule one can verify that r —>■ —2/k as r —)■ 3m. To 

analyse the behaviour of these curves as they approach the points Q and Q' let us consider r such 
that r = 3m + e . Then, one has that for small e > 0 that 

W{r) = ( —(— +^(r - 3m) + C>((r - Sm)^)) 

\r — 6in/ \m J 

where C\ and C 2 are numerical factors whose explicit form is not relevant for the subsequent 
discussion. Hence, to leading order lU(r) = Cjt^ where C is a constant depending on m only 
and p = 1/ \ph. Consequently, to leading order 

dr pCk€P~^ pCe 

^ ~ {MleP - kC) ~ MIcP-kC 

Therefore, since p < 2 one has that dr/de diverges as e —?► 0 so that the curves with 7 = 0 become 
tangent to the horizon as they approach Q or Q' —that is, they would have to become null to 
reach Q or Q'. This is analogous to the behaviour of the critical curve in the Schwarzschild 
spacetime pointed out in |23|, and the subextremal Reissner-Nordstrom spacetime in [42] — in 
contrast, in the extremal Reissner-Nordstrom spacetime one has ^ = 0 as e —>■ 0 as discussed in 

m ■ 


A. 4.2 Critical curves in the subextremal Schwarzschild-de Sitter spacetime 


For the subextremal case one could parametrise the roots of the depressed cubic (98) using Vieta’s 
formulae and choose some 7^1 for which there is at least one positive root. However, notice 
that fixing a value for 7 is equivalent to prescribe initial data for the congruence: 


t{f) = U, r(f) = r* 


dr 

df 




dt 

df 


7 

Fir*y 


( 100 ) 


Restricting our analysis to the static region < Vi, < Tc for which F(r*) > 0 and setting 

i = 0 . 

dr r* 

one gets 


7 = 


777 ), 


and condition (931 is equivalent to 


f (..) - F(r) = 

Sr 


where Q{r) is the polynomial 


Q{r) = r‘^ + r*r - 


6 m 

|A|r*' 


Notice that Q{r) can be factorised as 


50 
















Q(r) = (r - a_(r*))(r - a+(r*)), 


where 


In addition, observe that 


r* I , , L 24to \ 




Q{r*) 

> 0 

for 

r® < r* < rc, 

Qir*) 

< 0 , 

for 

rb <ri, < r@ 

Q{r*) 

= 0 , 

for 

r* = r@. 


\l/3 

where r 0 = j ■ In the extremal case one has r;, = Tc = r@ = 3m. The curve r = r@, as in 
the extremal case, will be called the critical curve. With the above notation the integral (92) can 
be then rewritten as 


T - = 


a 


(s - r*)(s - a_(r*))(s - a+(r*)) 


ds. 


( 101 ) 


To study the behaviour close to the critical curve consider r* = (1 + e)r@ For small e > 0 and 
considering s > r* one can expand the right hand side of equation (101) in Taylor series as 


T - r* = 


1 


3r|)Se^ 


+ 2 r@ys-r@ 2 (s - r@)3 


ds + C>(e^). 


Integrating we obtain 
2 


( 102 ) 


3r 


T - r* = - 


r + 2r, 


^arctanh + 2 In (^\/ r@(l + e) + ^r@(3 + e)^ - ^arctanh ^ 

+ 2 In ^\/r + + 2 r@^ — —+ 2 r 0 (l + 2 e) — —r@^'\/l + 2 r 0 -^- —\- 0{e^). 


(r® - ry 

As e —>■ 0 the last expression diverges —as is to be expected. The divergent term can be expanded 
for small e > 0 as 


arctanh ( 


'l±i =li„ 

3 + e/ 2 


-- + 4+ ^ + C>(e2) 
e o 


and the second term can be expanded as 

In ^@(1 + e) + y^r@(3 + e)^ = In (^(1 + V5)y/r^^ 
Hence, to leading order one has 


2^3 6^3 


+ 0(e3). 


^(0 = -^^n^ + f{r) + 0{e) 


where 


/(r) = f*+ 21 n ((1 + \/3),yr^)-^arctanh +2 In (v^+ i/r + 2 r@) -^rgi/l + 2 r@. 


V r + 2ri 


®' 
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Reparametrising respect to the unphysical proper time using (961 one gets 


2exp(^|A|/3/(r) + C)(e))eP 
1 - Kexp (\/|A|73/(r) + 0{e)) e? 

with p = 1/ -s/S. Thus one gets 

dr _ 2pexp (y/|A|/3/(r) + 0{e)) 

{l-exp{./\X\/3f{r) + 0{e))Kepf' 
Observe that since p < 1 then one has that dr/de diverges as e —>■ 0. 


B Appendix: The conformal evolution equations in the 
case K ^ 0 and reparametrisations 

In Section |4.2.2| we analysed the case k = 0 —this corresponds to a conformal boundary with 
vanishing extrinsic curvature. Nevertheless, as discussed in Section |2.4| k is a conformal gauge 
quantity arising from the conformal covariance of the conformal field equations. Consequently, 
it is of interest to analyse the behaviour of the core system in the case k ^ 0. For simplicity, 
in the remainder of this section, k will be assumed to be a constant on the initial hypersurface 
corresponding to r = 0. In first instance, we restrict our attention to |/t| > I and then discuss how 
to exploit the conformal covariance of the equations to extend these results for k G [—1,0) U (0,1]. 


B.l Analysis of the core system with k > 1 

We begin the discussion of this case observing that, for k > 1, one has that 0(t) > 0 and 
0(r) > 0 for r > 0. Using this simple observation and the core equations ( 73a| )-( [7^ we obtain 
the following: 


Lemma 14. For a solution to the core system (73a)-(73c I with initial data given by (74) and 
At > 1 one has that L(t) < 0 for r > 0. 

Proof. We proceed by contradiction. Assume that there exists 0 < < oo such that L{Tr) = 0. 

Without loss of generality we can assume that corresponds to the hrst zero of L{t). Since for 
At > 1 we have L(0) < 0 then by continuity it follows that L{ti) > 0 — L{tl) cannot be negative 
since this would imply that L(t) crossed the r-axis at some time t < tl but this is not possible 


since tl is the first zero of L{t). It follows then from equation (73c) that 


0 < L{tl) = -x{tl)L{tl) - ^ 0 ( rL ) 0 ( rL ). 

Since L{tl) = 0 and 0(tl) > 0, the last inequality implies that 4>{tl) < 0 but this is a contra¬ 
diction since we already know from Observation 1 that 0 (t) > 0 for any r. □ 


Observation 4. Using that 0(t) > 0 for At > I and r > 0 and that (/)(r) > 0 we obtain from 
equation (73c) the differential inequality 


L{t) < -x{t)L{t). 


Observing Lemma 14 we have that L[t) < 0. Thus, we can formally integrate the last differential 
inequality and obtain 

(103) 


L{t) < L{0) exp J x(s)ds^ . 


We now show that the function x(t) which is initially positive must necessarily have a zero. 
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Lemma 15. For a solution to the core system (73a)-(73c) with initial data given by (74) and 
K > 1 there exist 0 < < oo such that x(t^) = 0. 

Proof. We proceed again by contradiction. Assume that x(t) never vanishes. Since x(0) = k > 0 
then x(r) > 0 for r > 0. From Lemma 14 we know that L{t) < 0. In addition, we know that 
0(t)(/)(t) > 0. With these observations equation ( |73b ) gives 

x { t ) < -X^(t) for T > 0. 

Since we are assuming that x(t) never vanishes then 

X{r) 


X^{t) 


< - 1 . 


Integrating from 0 to r > 0 and using the initial data (74) we get 

I 


x{t) < 


T -\-l/h 


for r > 0. 


(104) 


In a similar way, we can consider equation (73b[) and obtain the differential inequality 


X('r) < for r > 0. 


Using now equation (751 we get 


X < —TO0(r) exp ^—3 J x(s)ds^ for r > 0. 
Integrating the from 0 to r > 0 we get 

x(r) < K — m J 0(s) exp ^—3 J x(s')ds'^ ds for r > 0. 


(105) 


On the other hand, integrating expression (1041 we have 


Consequently, 


Integrating we get 


/ x(s)ds < In {kt + 1). 
Jo 

-me(r)exp(-3^ 


e(s)exp (-3^ x(3')*') * < “ 0 ' 


Substituting the above result into the inequality (105) we obtain 

+ ln(KT + 1) — 1^ . 


m \X\ 

x{r) <n-ixx2\hr 


1 


V 3 \ {kt + 1) 

The right hand side of the last expression becomes negative for some sufficiently large r. This is 
a contradiction as we have assumed that x(r) never vanishes and x(0) >0. □ 
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Observation 5. Combining Lemma 14 and Observation 1, we conclude that L(t) < 0 and 
0(r)^(T) > 0 for T > 0. Using these properties in equation (73b) we get 

x(t) < 0 for T > 0. 


Thus, x(t) is always decreasing. From Lemma 15 we know that there exists a finite > 0 such 
that xi^x) ~ 0. Then, by continuity, for any t > we have that xi''') < 0 • 

With this last observation we are in the position of proving the main result of this section: 


Proposition 1. There exists 0 < < oo such that the solution of (73a)-(73c I with initial data 

given by (74) and k> 1 satisfies 


X 


n- 


Proof. From Lemma 15 we know there exists a finite for which x('r) vanishes. By Observation 
5, we have that x('''o) < 0 for any Let xo = x(to) < 0. We can assume that xo is finite, 

otherwise there is nothing to prove. Now, using Lemma [l4] and that 0(r)(^(T) > 0 we get 


X('r) < -X^(t) for T > 0. 


Since we know that x(''‘) < 0 for any t > then 


x{t) 

X^{t) 


< - 1 . 


Integrating form r = rq to t > to we get 


(106) 


From inequality (106) we can conclude that x('^) — oo for some finite time < rq — 1/xo- 

Additionally, observe that rq — 1/xo > "rq > 0 since XO < 0 • Now, given that x —oo as 
T —> T^ it follows from equation (75) that <()—)■ oo as r —)■ . Similarly, from inequality ( [Tost 
and that L(0) < 0 it follows that L —)■ —oo as r —>■ . □ 


Remark 10. A plot of the numerical evaluation of the solutions to the core system (73a)-(73c) 
with initial data (74) in the case k > 1 can be seen in Figure]^ 


B.2 Analysis of the core system with k < —1 

In this section we use a similar approach to that followed in Section p3.1| to show that the fields in 
the core system diverge for some finite time if k < — 1. An interesting feature of this case is that, 
assuming one knows that there exists a singularity in the development, there exists an a priori 
upper bound for the time of its appearance —namely, the location of second component of the 
conformal boundary at r = 2/|«:|. As a byproduct of the analysis of this section an improvement 
of this basic bound is obtained. 


An important remark concerning the case k < —1 is that if r G [0, 1/|k|] then both 0(t) and 
0(t) are non-negative. Based on this observation our first result is: 


Lemma 16. If k < — 

satisfies L(t) <0 for t 


1 then the solution to the core system (73a)-(73c) with initial data (74) 

e [0,1/|«:|]. 
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Figure 9: Numerical solutions of the core system (73a)-(73c) with initial data given by (74) in 
the case k = 2 and |A| = 3, m = l/3-\/3- The solid line describes the evolution of </), the dashed 
line that of x s-nd the dotted line that of L. One can observe the formation of a singularity at 
T « 2.6392. 


Proof. We proceed by contradiction. Assume that there exists 0 < < l/|At| such that L{ti,) = 

0. Without lost of generality we can assume that tl is the first zero of L(t). Since T(0) < 0 for 
K < —1 then by continuity L{Tif) > 0. Therefore, proceeding as in Lemma 14 one gets from (73c) 

0 < T(tl) = -x{tl)L{tl) - ^0(rL)0(rL) for r e [0, 1 /|k|]. 


Since L{tl) = 0 and ©(r^) > 0 the last inequality implies that < 0. This is a contradiction 

since (/)(r) > 0 —cfr. Observation 1. □ 


Lemma 17. If k < —1 then the solution to the core system (73a)-(73c) with initial data (74) 
satisfies xi^) < 0 for r G [0, 1 /|k|]. 

Proof. Again, we proceed by contradiction. Assume that there exists 0 < < 1 /|k| such that 

xi^x) = 0. Without lost of generality we can assume that is the first zero of x(t). Then, by 
continuity, we have that x(t^) > 0. Using equation (73b) one has 


0 < x(w) = -X(W) + L{t-x) - 


for r G [0, I/\k\]. 


Therefore, since x(w) = ^ one has 


L{t^) > ^Q{Tx)(t>{Tx) > 0 . 


This is a contradiction since by Lemma 16 we know that Lir) < 0 for r G [0, 1 /|k|]. 
Observation 6. Proceeding as in Observation 4 one readily has that for k < — 1 
T(r) < L(0) exp y x(s)ds^ for rG(0,l/|K|]. 

This last observation is used, in turn, to prove the main result of this section: 


□ 
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Figure 10: Numerical solution of the core system (73a)-(73c) with |A| = 3, to = l/3-\/3 in the 
case K = —2. The solid line corresponds to </>, the dashed line to % and the dotted line to L. One 
can observe a singularity at t « 0.4203. 


Proposition 2. If k < —1, then for the solution of (73a)-(73c) with initial data (74) there exists 
0 < < \/\k\ sueh that 

x(t) —>■ —00, L{t) —)■ —00, and (fir) —>■ oo as t ^ 


Proof. Consider equation (73b) on the interval r S [0, 1 /|k|]. Using Lemma 16 we know that 
L(t) < 0. This observation and the fact that > 0 leads to the differential inequality 

Xir) < for T e [0,1/|«:|]. 


Since by Lemma 17 we know that x('^) 7^ 0 ^ [0? 1/1^1] rewrite the last expression 

as 

X{r) 




< —1 for T G [0, 1/\k\]. 


Integrating from t = 0 to 1/\k\ and using the initial data (74) we get 

x{t) 


(107) 


From inequality ( |107 ) one concludes that x(t) — >■ — oo for some 0 < < 1/\k\. Finally, using 

Observation 6 and Observation 1 one concludes that L{t) —>■ —oo and ^(t) —>■ oo as t —?> for 
some 0 < < l/|/c|. □ 

Notice that this upper bound for the location of the singularity is not trivial and improves the 
basic bound r < 2 /|k| given by the location of the second component of the conformal boundary. 


Remark 11. A plot of the numerical evaluation of the solutions to the core system (73a)-(73c) 
with initial data (74) in the case k < —1 can be seen in Figure [T^ 


B.3 Exploiting the conformal gauge 

In Lemma |12| we have shown that if = 0 then the evolution equations imply, in particular, 
fx = 0. Due to the spherical symmetry Ansatz, the component fx is the only potentially non-zero 
component of /. Thus, one concludes that f — 0. In Section [B.3.1 we will exploit this feature 
of the Weyl connection to extract further information about k and s. These results are used in 
Section [B. 3. 2| to discuss the conformal gauge freedom of the extended conformal field equations 
and the role played by reparametrisations of conformal geodesics. 
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B.3.1 The relation between the Weyl and Levi-Civita connections 


As discussed in Section 2.2.1 the Weyl connection V expressing the extended conformal field 
equations is related to the Levi-Civita connection V of the unphysical metric g via the 1-form f. 
If / vanishes then V = V. Exploiting this simple observation we obtain the following results: 


Lemma 18. // / = 0 then the conformal gauge conditions (12) and (13) imply that s = 0. 
Moreover, s is constant along the conformal geodesics. 


Proof. As discussed in Section |2.2.l| if / = 0 then Lab = Lab and fa^b = Ta'^b- Using the 
conformal gauge condition (12) it follows that Loa = 0 and Lo^b = 0. Now, the standard 
conformal held equations (124d) and (124e) in Appendix [P] give 


VqVq© -|- ©Lqo ~ srjoo — 0, 
Vos = -LobV^©. 


(108a) 

(108b) 


Using Loa — 0 and ro“b = 0 in equation (108a) one concludes © = s. Similarly, from equation 
(108b) one gets s = 0. Therefore s is constant along the conformal geodesics. □ 


Remark 12. In the asymptotic initial value problem the initial value of s is given by s* = 
—see equation (34a). Thus, if f = 0 then s = ■\/|A|/3k along the conformal geodesics. 


Finally, one has the following: 

Lemma 19. In the asymptotic initial value problem, if f = 0, then the conformal gauge conditions 


(12) and (13) together with the conformal Einstein field equations imply that ei^n) — 0 — that is, 
K is a constant. 


Proof. Using / = 0 and the gauge conditions (12) we get from the conformal held equation (124g) 
that 

60s - 30^ -f 3(5*^ei©ej© = A. (109) 


Using Lemma 18 we have s = 0. Therefore, substituting 0(r) = 0 *t(1 -|- kt/ 2) into equation 
(109) and recalling 0* = a/|A|/ 3 we obtain 


ei{K)ej{K) = 0 . 


Observe that the last equation is trivially satished on as r = 0. Oh the initial hypersurface, 
where t ^ 0, the last equation implies 


6^^ ei{K)ej{K) = 0 . 


Therefore, we conclude that ei{K) = 0. 


□ 


B.3.2 Changing the conformal gauge 


The analysis of the core system given in Sections |B.1[ m and Section |4.2.2| covers the cases 
for which |k| > 1 and k = 0. As a consequence of the conformal covariance of the extended 
conformal Einstein held equations one has the freedom of performing conformal rescalings and 
of reparametrising the conformal geodesics —thus, ehectively changing the representative of the 
conformal class [g\ one is working with. This conformal freedom can be exploited to extend the 


analysis given in Sections B.l and B.2 to the case where k G [—1, 0) U (0,1]. 

Following the discussion in the previous paragraph, any two spacetimes {M,g) and {M,g) 
with g — Q^g and g = Q^g representing two solutions to the extended conformal Einstein held 
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equations for different choices of parameter k are conformally related. From Lemmas and we 
have that 


e{r) = ^r[l+l^r), 0(f) = ^f( 


1 + 


with 


( 110 ) 


( 111 ) 


CT + d 

The free parameter b in the fractional transformation of Lemma has been set to 6 = 0 in 
order to ensure that 0 and 0 vanish at r = 0 and f = 0, respectively. Thus, the conformal 
boundary ^ is equivalently represented by the hypersurfaces with r = 0 or f = 0. As g and g 
are conformally related one can write 


— cici-i 


g = oj'^g with oj = 00 


Using the relations in (110) and (111) we obtain, after a calculation, that 

gRt 


}{t) = 


2(ct + d) 


a ( 1 + 


(cT + d) ( 1 + -KT 


( 112 ) 


The conformal transformation law for the field s can be seen to be given by 

S = UJ~^S + LJ~^VcU!V^0 + ^lj~^0VcGjV^uj. 

As discussed in Section |B.3.1[ in the analysis of the extrema l Sc hwa rzsc hild-de Sitter spacetime 


one can assume that d^pK = 0 and / = 0. Now, Propositions 


18 


and 


19 


imply that s = ^J\\\/?>K 


and s = are constant. Exploiting this observation, the transformation law for s can be 

read as an equation for w —namely 


0w^ + 2w0w + uj^s — uj^s = 0. 


Substituting expression (112) into equation (113) one gets the condition 

2c + aR — dn = 0. 


(113) 

(114) 


One can read equation (114) as the transformation law for k so that 

dn — 2c 


In order to have a meaningful transformation law between r and t, neither a nor d can vanish. 


Substituting equation (114) into the reparametrisation formula (111) and expression (112) one 

:a;*. T 

w(t) = 


can observe that a/d actually corresponds to w(0) = w*. Therefore, one has that 

2w*t , , 4a;* 


-(t) = 


(a;*«; — k)t — 2 ’ 


((a;*K - k)t - 2) 


2 ■ 


(115) 


From the last expression one can identify w* = d;(0) = uji,{u)i,R — k). In addition, notice that 
f —>■ oo and a; —)■ oo as r —>■ 2/(a;*lt — k). Therefore, the hypersurface defined by r = 2/(a;*K — k) 
is at an infinite distance from the conformal boundary as measured with respect to the g-proper 
time. 


Remark 13. An alternative approach to deduce equations (114) and (115) is to write Q(t{t)) = 
and use equations (110) and (111) to identify k and a;. 
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C Appendix: Cartan’s structure equations and space spinor 
formalism 

In this appendix we give a brief discussion of Cartan’s structure equations and the space spinor 
formalism. 

C.l Cartan’s structure equations in frame formalism 

Consider a ft,-orthonormal frame {e*} with corresponding coframe {u;*}. By construction, one 
has (a;*,ej) = Si^. The connection coefficients of the Levi-Civita connection D oi h respect to 
this frame are defined as 

{u)^,Diek) =li^k- 

As a consequence of the metricity of D it follows that jijk = —"fikj- The connection form is 
accordingly defined as 

= 7i^fc A uj\ 

With these definitions, the first and second Cartan’s structure equations are, respectively, given 
by 


du;® = —A , (116a) 

d7%- = -7'feA7^- + 0%-, (116b) 


where r2®j is the curvature 2-form defined as 

Aw*. 


C.2 Basic spiuors 

In the space spinor formalism, given a spin basis where a=o .17 any of the spinorial fields 

appearing in the extended conformal Einstein field equations can be decomposed in terms of basic 
irreducible spinors. The basic valence-2 symmetric spinors are: 


xab = 

The basic valence 4 spinors are given by 


_ 1 11 
Vab = ’ 


SACXbD + eBDXAC, 


eAcVBD + eBnVAC, 


^ ^ ^ 0^ 0 
^ab = -^Ha es ) • 


iACZBD + SBDZAC, 


hABCD = —^A(C^D)B, ABCD = ^(A^^ ^B^ 


(117) 


(118a) 

(118b) 


In the last expression iABCD)i in^jeates that an i number of indices are set equal to 1 after 
symmetrisation. Any valence 4 spinor Cabcd with the symmetries C{ab){cd) can be expanded 
in terms of these basic spinors. One has the identities 


2 ^4 ^0 

X(abXcd) = 2 e ABcn, y(ABVcD) = 2 ^ abcb, Z(^abZcl>) = 2 ^ abcb, 


1 


VABXcD = —e ABCD 


2i/2 
1 


{iACUBD + ^BDUAc), 


ZABXcD = tABCn + —J={£ACZBD + ^bdzac) 


1 2 

UABZcD ~ ~ 2 ^ ABCD 


2i/2^ 

1 . ^ Iz. 

-—f=\^ACXBD + ^BDXac) — -^hABCD- 
4 V 2 6 


(119a) 

(119b) 

(119c) 

(119d) 
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Another set of identities used in the main text is given by 


xabx 


AB 


= 1, xabU^^ = 0, xabz^^ = 0, zabz^^ = 0, habz^^ = (120a) 


xa^xbq = 2^ab, Ua^xbq = -^Uab, za^xbq = -^zab, Va^Vbq = 0,(120b) 


Va'^zbq = 


1 


2^/2 


XAB 


^CAB, Za^ZbQ 


yCD _ 


1 


1 


= -^xab, 


^abcdv'"^ — qVab, 


’^ABCD 

These identities and a more exhaustive list has been given in m- 


2 OD _ ^ , 

^ABCD^ — q^AB- 


(120c) 

(120d) 


C.3 Cartan’s structure equations in spinor form 

Let denote a Hermitian spinor with normalisation taa' = 2. Consider an adapted 
spin dyad such that the matrix representation of is given by the identity 2x2 matrix. 

The spatial Infeld-van de Waerden symbols are related to the usual Infeld-van der Waerden via 

( 121 ) 


..r i — ^ B' 

(XAB = T(b (XA)B' 


Equivalently, one has 


^AB _ (B „A)B' 

(J 2 — ^ B'^ i' 

The matrix representation of the spatial Infeld-van der Waerden symbols is given by 


(o‘ 1 


AB _ 1 /^—1 0 

" ^ = 7 iU 1 


2 _ 1 fi 0 

=75Vo i 

AB _ J_ f-i 0 

"“y2V0 i 


^ 3_ WO A 

=7ill Oj’ 

AB _ fO 1^ 

" “ 72 U oj ■ 


Thus, the space spinor counterpart of coframe and connection coefficients can be obtained suc¬ 
cinctly by contraction with the spatial Infeld-van der Waerden symbols as and 

')ab'"^EF = k<x''ABfXj'^^cr^EF- With these definitions the spinorial version of the Cartan 
structure equations is given by 


where 


= —j^B A o)'®'® — 7 '^b a 
d7^B = -7 ^e a 7^b + il^B, 


^,A — ^ aq , nn 

7 B = bq<^ , 


(122a) 

(122b) 


and is the spinor version of the curvature 2-form, with 


n^B = 


bcdef(^^^ a . 


In the last expression the spinor rABCDEF can be decomposed as 


XABCDEF = ( ^SaBCD — -^xhABCE ) ^DF 


-j^SABDF — J^rhABDF ) ^CE 


where sabcd and r correspond to the space spinor version of the trace-free part of the Ricci 
tensor and Ricci scalar of h, respectively. 

To relate the previous discussion with the basic spinors xab, Uab and zab, observe that 
using (|117|) and (|121|) one obtains that 


(XAB =—ZAB —UAB, 


(XAB^ = K^AB — VAb), (XAB^ = XaB, 


^ab 


2 = K — ZAB + Vab), 


= -x^B 


(123a) 

(123b) 
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D Appendix: The frame conformal Einstein field equations 


The tensorial (frame) version of the standard vacuum conformal Einstein field equations are given 
by the following system —see e.g. [iniiniiiaiTS]: 


(124a) 

(124b) 

(124c) 

(124d) 

(124e) 

(124f) 

(124g) 


= 0 , 

Ved'^ab/ = 0, 

V c^db V dLbc V bed — 0 
V4" ^d/ab ^9ab — O 5 
VaS + TaeV'^.n. = 0, 

R abd P abd — 0; 

65s - 3VaSV“5 = A, 


where Sa'^b is the torsion tensor, given in terms of the connection coefficients, as 

b^c = ^b] (d'a b a)^ci 

Lab is the Schouten tensor; 5 is the conformal factor and s is a concomitant of the conformal 
factor defined by 

In addition, p“bcd is the algebraic curvature and R'^dab is the geometric curvature. 


P“bcd = ^d'^bed + 2(5“[eTd]b - 5b[c'^d]“), 

R^ dab = ea(r5‘^d) — eb(ra‘^d) + r/‘^d(rb'^a — ha-^b) 
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